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ABSTRACT 


A  new  model  is  reported  for  the  prediction  of  the 
electrical  conductivity  of  polymeric  materials  filled  with 
conducting  particles,  short  fibers,  and/or  platelets.  The 
model  takes  into  account  chain  and  network  formation  by 
contacting  filler  units  which  provide  continuous  conduction 
paths.  The  incorporation  of  this  important  structural 
feature  shows  that  sharp  transitions  in  conductivity  occur 
within  a  narrow  range  of  filler  concentrations.  The 
critical  concentrations  associated  with  these  transitions 
establish  the  minimum  concentration  of  conducting  filler 
required  to  achieve  conducting  behavior. 

Preliminary  evaluations  of  the  model  indicate  that 
it  can  be  used  to  rank  candidate  materials.  Predictions  of 
absolute  levels  of  conductivity  require  knowledge  of  (i) 
experimentally  determined  values  for  the  conductivity  of 
each  of  the  components,  (ii)  the  volume  fraction  (or  weight 
fraction)  of  each  of  the  components,  (iii)  an  independent 
characterization  of  the  average  state  of  orientation  of  the 
conducting  components,  and  (iv)  an  independent  characteriza¬ 
tion  of  the  average  aspect  ratio  associated  with  the  filler 
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shape.  The  characteristic  parameters  of  the  contact 
probabilities  (determined  by  the  orientation  and  aspect 
ratio  distributions)  are  particularly  sensitive  to  pro¬ 
cessing  conditions.  Consequently,  variations  in  processing 
procedures  can  have  a  pronounced  effect  on  the  electrical 
properties.  This  result  emphasizes  the  importance  of 
conducting  a  microstructural  characterization  program  in 
concert  with  electrical  characterizations.  A  procedure  is 
identified  for  determining  the  essential  microstructural 
data . 


The  routine  application  of  the  model  relationships 
as  a  screening  tool  for  candidate  material  systems  requires 
tedious  com.putations .  An  interactive  computer  program  was 
developed  to  conduct  these  routine  calculations.  A  user's 
guide  with  illustrative  examples  is  given  in  the  appendix. 
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INTRODUCTION 


The  rapidly  evolving  electronics  industry  has 
prompted  the  need  for  light  weight,  low  cost  materials  that 
can  be  formed  by  existing  processing  procedures  into 
housings  that  possess  electromagnetic  shielding  charac¬ 
teristics  resembling  metals.  Considerable  interest  has 
developed  in  the  use  of  corrosion  resistant  polymeric 
materials  filled  with  conducting  particles,  fibers,  and/or 
platelets.  The  word  "metalloplastic"  [Davenport,  1981]  has 
been  coined  to  distinguish  this  class  of  heterogeneous 
materials  from  the  emerging  class  of  homgeneous  semi¬ 
conducting  and  conducting  polymeric  materials. 

One  of  the  key  material  properties  required  of 
metalloplastics  in  these  applications  is  electrical  conduc¬ 
tivity  approaching  the  range  of  metals.  The  wide  range  of 
possible  compositions  and  filler  geometries  (e.g.,  fibers, 
particles,  platelets)  precludes  a  total  reliance  on  direct 
laboratory  characterization  of  the  directionally  dependent 
electrical  properties  of  these  materials.  Accordingly,  an 
objective  of  this  work  was  to  develop  models  to  predict  the 
electrical  properties  of  metalloplastics  with  the  intent  of 
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providing  a  screening  tool  for  the  identification  of 
composite  systems  which  are  conducting  and  exhbit 
acceptable  properties . 

It  is  evident  that  such  a  model  will  contain 
parameters  associated  with: 

o  the  electrical  properties  of  the  components; 

o  the  composition  (e.g.,  volume  fraction  or 
weight  fraction) ; 

o  the  shape  of  the  conducting  filler;  and 

o  the  relative  orientation  and/or  packing 
geometries  of  the  conducting  fillers. 

It  can  be  anticipated  that  the  computational  complexities 
introduced  by  the  necessity  of  treating  such  an  extended  sat 
of  parameters  will  require  computer  programs  for  the  imple¬ 
mentation  of  the  model.  Accordingly,  a  second  objective  of 
this  work  was  to  develop  interactive  computer  programs  to 
facilitate  the  use  of  the  models  as  a  screening  tool. 


CLASSIFICATION  OF  TRANSPORT  PROPERTIES 


Electrical  conductivity  belongs  to  the  general  class 
of  transport  properties  which  include; 

o  thermal  conductivity 
o  dielectric  constants 
o  magnetic  permeability 
o  diffusivity 

As  noted  by  \7u  and  McCullough  [1977],  among  others,  each  of 
these  properties  is  governed  by  linear  constitutive  rela¬ 
tionships  which  connect  a  response  field  to  an  applied  field 
through  an  appropriate  descriptor  of  the  material.  These 
descriptors  (e.g.,  electrical  and  thermal  conductivity)  are 
the  key  "intrinsic"  material  parameters  that  are  required 
for  the  evaluation  of  the  performance  of  various  designs  of 
fabricated  parts  of  complex  shapes.  In  particular, 
computer-aided  design  techniques  aimed  at  achieving  such 
evaluations  require  these  material  parameters  as  input.  Such 
an  evaluation  would  include  the  processing  of  position 
variant  fiber  distribution  and  orientation  information 
[Quigley,  1982]  into  a  format  for  a  finite  element  analysis 
package  such  as  ADINAT  [1982]. 
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Although  the  underlying  physical  processes  differ 
for  each  of  the  transport  properties,  the  mathematical 
structure  for  model  relationships  that  connect  the  proper¬ 
ties  and  microstructure  of  the  components  of  these  composite 
systems  to  the  effective  (average)  bulk  descriptors  will  be 
similar.  Consequently,  it  can  be  anticipated  that  models 
developed  for  electrical  conductivity  will  find  application 
in  predicting  other  transport  properties  (e.g.,  thermal 
conductivity  and  dielectric  constants) . 

The  general  structure  of  linear  transport  relation¬ 
ships  is  of  the  matrix  form: 


where  the  G's  are  components  of  the  applied  field  (e.g., 
electrical  fields,  temperature  gradients,  chemical 
potentials)  and  the  F's  are  components  of  the  corresponding 
response  field.  The  Kj[j's  are  the  appropriate  descriptors 
of  the  material  (e.g.,  thermal  conductivity,  electrical  con¬ 
ductivity,  and  diffusivity,  respectively) .  Symmetry 
requires  that  Ki j  =  Kji  .  These  relationships  can  be 
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simplified  by  restricting  attention  to  the  important  class 
of  composite  materials  which  exhibit  orthotropic  symmetry 
characterized  by  three  mutually  perpendicular  planes  of 
material  symmetry.  If  the  "1",  "2",  and  "3"  directions  are 
selected  as  the  principle  axes,  only  the  diagonal  terms  of 
[K]  survive,  so  that  Equation  (1)  reduces  to 


K 
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0 


0  K22 


0  0 


(2) 


where  Kn,  K22»  K33  can  be  selected  to  refer  to  the 

"longitudinal",  "transverse"  and  " through-the-thickness " 
properties  of  a  rectangular  panel  specimen.  The  behavior  of 
a  transversely  isotropic  material  can  be  obtained  as  a 
special  case  by  setting  any  two  of  the  K's  equal  (e.g., 

K22  =  K33);  isotropic  material  symmetry  is  obtained  by 
equating  all  K's,  viz.,  Kn  =  K22  =  K33  .  Behavior  in 
directions  other  than  the  principle  axes  can  be  obtained  by 
projecting  the  K's  through  well  known  rotational  transfor¬ 
mations.  As  a  consequence  of  this  reduction,  the  double 
index  "ji"  can  be  replaced  by  a  single  index,  i.e., 

Kji  =  Kj  . 
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Based  on  this  general  formulation,  models  developed 
for  the  various  K’s  (representing  electrical  conductivity) 
may  be  extended  to  the  prediction  of  other  transport  proper¬ 
ties  by  replacing  the  appropriate  modification  of  the  K's 
with  the  transport  property  of  interest. 


REVIEV7  OF  TRANSPORT  MODELS 


The  literature  is  replete  with  relationships 
proposed  to  predict  various  transport  properties.  Models 
for  thermal  conductivity  are  reviewed  by  Progelhof  et  al 
[1976];  models  for  electrical  and  magnetic  properties  are 
reviewed  by  Hale  [1976].  As  with  thermomechanical  proper¬ 
ties  [McCullough,  1971],  the*various  approaches  to  model 
development  may  be  broadly  classified  as  "Mechanics  of 
Materials,"  "Self-Consistent  Field,"  and  "Bounding." 


The  mechanics  of  materials  approach  directs 
attention  to  simple  geometries  and  introduces  gross  approxi¬ 
mations  to  the  response  fields.  For  example,  the  simplest 
of  such  models  would  correspond  to  continuous  conducting 
fibers  all  aligned  parallel  to  the  "1"  direction.  The  con¬ 
ductivity  in  the  "1"  direction  is  obtained  by  considering 
that  the  conducting  fibers  are  connected  in  parallel  with  an 
"insulating"  matrix  so  that 


Ki  =  v^K^  t 


V  K 
m  m 


(3a) 


where  the  subscripts  denote  the  fiber  "f"  and  matrix  "m" , 
respectively;  Vf  is  the  volume  fraction  of  the  fiber  with 
=  1  -  Vf  . 
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The  conductivity  in  the  "2"  (or  "3")  direction  is 
obtained  by  considering  that,  in  the  transverse  direction, 
the  conducting  fibers  and  insulating  matrix  respond  through 
a  series  connection  so  that 


(l/Kp  =  (Vj/Kj)  +  (V_^/K^) 


or 


P-,=  v  p  +v  p 
2  f  f  m  m 


(3b) 


The  "rule  of  mixtures"  expressed  by  Equations  (3)  has  been 
verified  as  an  appropriate  model  for  predicting  the  longi¬ 
tudinal  transport  properties  for  the  special  case  of  con¬ 
tinuous  and  aligned  fibers  [Equation  (3a)].  However,  these 
simple  models  fail  to  predict  the  transverse  properties 
[Equation  (3b)],  and  are  inadequate  for  particulate  and 
short  fiber  filler  geometries. 

In  contrast,  the  self-consistent  field  approach 
introduces  gross  approximations  in  phase  geometry  and 
simplifies  the  representation  of  the  response  field  in  an 
attempt  to  obtain  rigorous  solutions  to  boundary  value 
problems.  These  models  generally  overestimate  transport 
properties.  The  bounding  approach  avoids  the  difficulties 
introduced  by  complex  microstructures  by  directing  attention 
to  establishing  upper  and  lower  bounds  on  the  transport 
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properties  rather  than  predicting  the  properties.  The 
results  from  these  bounding  approaches  can  serve  as  a  prac¬ 
tical  guide  to  material  behavior  only  if  the  upper  and  lower 
bounds  are  sufficiently  close  so  as  to  bracket  the  actual 
behavior  within  experimental  error.  In  general,  the 
transport  properties  of  filler  and  polymer  differ  by  several 
orders  of  magnitude  so  that  the  upper  and  lower  bounds  are 
too  far  apart  to  provide  useful  results. 

All  of  these  approaches  have  their  counterpart  in 
developing  models  to  predict  the  elastic  constants  of  com¬ 
posite  materials.  Indeed,  several  workers  have  noted  the 
similarity  between  the  structure  of  Equation  (2)  to  the 
shear  modulus  "block"  of  the  elastic  constant  array  and  have 
proposed  that  transport  properties  can  be  predicted  by 
replacing  the  shear  modulus  in  various  elastic  models  by  the 
transport  properties  [e.g.,  Aston  et  al,  1969].  This 
analogy  appears  to  give  reasonable  results  for  continuous 
fiber  composites;  it  fails  to  describe  the  behavior  of  short 
fiber  and  particulate  filled  materials. 

While  some  models  appear  adequate  for  predicting 
certain  transport  properties  for  restricted  filler 
geometries,  no  model  provides  the  generality  demanded  for 
the  prediction  of  the  entire  class  of  transport  properties 
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for  a  variety  of  filler  geometries.  As  pointed  oat  by 
Nielsen  [1974],  among  others,  the  deficiencies  of  the 
various  models  stem  from  the  neglect  of  important  structural 
features.  Noting  that  the  Halpin-Tsai  relationship  [Aston 
et  al,  1969]  was  a  convenient  form  for  representing  upper 
and  lower  bound  relationships,  as  well  as  results  from  self- 
consistent  field  models,  Lewis  and  Nielsen  [1970]  introduced 
modifications  in  the  Halpin-Tsai  relationship  to  partially 
account  for  the  shape  and  packing  geometry  of  inclusions. 
These  workers  argue  that  the  volume  fraction  of  filler  can¬ 
not  exceed  a  value  associated  with  close  packed  geometries 
(as  determined  by  the  shape  of  the  filler).  Accordingly, 
the  volume  fraction  of  filler  was  adjusted  to  account  for 
the  limitation  on  the  maximum  volume  fraction.  Estimates  of 
the  modification  parameters  for  various  collections  of 
spheres  and  rods  are  given  by  these  workers.  Predicted 
values  for  elastic  moduli  and  thermal  conductivity  obtained 
from  this  modification  are  in  better  agreement  with  observed 
behavior  than  predictions  obtained  from  models  which  ignore 
the  influence  of  the  microstructure.  McGee  and  McCullough 
[1981]  have  shown  that  the  Lewis-Nielsen  modification  of  the 
Halpin-Tsai  relationship  tends  to  overestimate  elastic 
moduli  at  high  filler  concentration.  These  workers  report  a 
new  model  that  gives  good  agreement  between  predicted  and 
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Observed  elastic  moduli  for  particulate  filled  systems  as 
well  as  porous  materials.  Bradbury  and  Bigg  [1980]  show 
that  while  the  thermal  conductivity  of  aluminum  and  copper 
I  filled  epoxies  are  reasonably  approximated  by  the  Lewis- 

Nielsen  relationship,  the  electrical  conductivity  exhibits  a 
sharp  transition  that  is  not  accounted  for  by  the  Lewis- 
Nielsen  equation,  or  any  other  available  relationship. 

Gurland  [1966]  proposed  that  the  observed  sharp 
transitions  in  electrical  conductivity  of  metal  filled 
polymers  was  a  consequence  of  the  formation  of  chains  of 
contacting  fillers  which  provide  continuous  paths  for  con¬ 
duction.  The  relative  frequency  of  such  continuous  paths 
should  be  dependent  upon  the  shape  of  the  filler  and  the 
volume  fraction  of  filler.  The  Lewis-Nielsen  relationship 
requires  that  such  chain  formation  occur  only  in  the 
vicinity  of  the  maximum  packing  fraction  of  filler. 

Evidently,  significant  chain  formation  can  occur  at  con¬ 
siderably  lower  volume  fractions  of  filler.  This  effect  is 
amplified  when  tne  transport  phenomena  can  be  distinguished 
as  volume  (e.g.,  thermal  conductivity)  versus  surface  (e.g., 
electrical  conductivity)  transport  processes. 


FORMULATION  OF  NETWORK  MODEL 


In  view  of  Gurland's  [1966]  and  other  similar 
proposals,  it  is  reasonable  to  assume  that  chain  and  network 
formation  is  the  origin  of  the  sharp  transitions  observed  in 
the  electrical  conductivity.  The  qualitative  significance 
of  chain  and  network  formation  is  illustrated  in  Figure  1 
for  a  particulate  filled  polymer.  The  resistivity  of  the 
particulate  filler  is  taken  as  =  10~6  ohm-cm 
[Kf  =  10^  (ohm-cm) ~1];  the  polymeric  binder  has  a  resistiv¬ 
ity  of  =  10^^  ohm-cm  [Kj^  =  10“1®  (ohm-cm)“l].  As  illus¬ 

trated  in  Figure  1 ,  the  resistivity  of  the  composite  may 
range  by  several  orders  of  magnitude  from  that  of  an 
insulator  (p  =  10+14  _  lo+l®  ohm-cm)  to  that  of  a  conductor 
(p  =  10~2  _  io~3  ohm-cm)  as  the  volume  fraction  of  the  con¬ 
ducting  filler,  Vf  is  increased.  The  transition  from  insu¬ 
lator  to  conductor  can  be  rationalized  by  the  following 
structural  arguments. 

Structural  Basis 

At  low  volume  fractions,  the  conducting  fillers  are 
isolated  by  the  insulating  matrix  phase.  Consequently,  the 
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Figure  1.  Illustration  of  the  role  of  chain  and  network 

formation  in  establishing  the  electrical  reisis- 
tivity,  p  ,  of  a  metal  particle  filled  polymer. 
The  open  circles  represent  isolated  metal 
particles.  The  closed  circles  represent  contact 
ing  particles  that  form  chains  and  networks  of 
continuous  conduction  paths  as  the  volume 
fraction,  Vf  ,  of  the  particles  is  increased. 
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composite  tends  to  exhibit  the  high  resistivities  of  the 
matrix  (p  =  ohm-cm)  and  behaves  as  an  insulator.  In 

this  insulating  plateau  region,  the  resistivity  slightly 
decreases  with  increasing  concentration  of  conducting 
particles  as  approximated  by  Equation  (3b).  As  the  con¬ 
centration  of  conducting  particles  is  further  increased,  the 
probability  of  particle-particle  contact  is  enhanced.  These 
contacts  promote  the  formation  of  bead-lihe  chains  which 
provide  continuous  conduction  paths  that  mimic  the  behavior 
of  filaments  of  the  the  conducting  materials.  The  transi¬ 
tion  from  bead— like  to  fiber-like  behavior  occurs  at  a 
critical  volume  fraction,  Vf  =  ,  for  which  the  average 

number  of  particle-particle  contacts  begins  to  exceed  "1". 
Based  on  Gurland's.  observations,  this  critical  concentration 
is  in  the  vicinity  of  Vf  =  0.3.  Further  increases  in  the 
volume  fraction  concentration  of  the  conducting  particles 
causes  a  progressive  increase  in  the  average  chain  length  of 
contacting  particles  so  that  the  behavior  rapidly  tends 
toward  the  limiting  behavior  of  a  continuous  fiber  system. 
The  structure  associated  with  the  region  intermediate 
between  insulating  and  conducting  behavior  is  illustrated  in 
Figure  1  by  the  closed  circles  which  represent  contacting 
particles  participating  in  chains  of  varying  conduction  path 
lengths . 
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A  continued  increase  in  the  concentration  of 
conducting  particles  causes  an  increase  in  the  average 
number  of  contacts.  Eventually,  the  average  number  of  con¬ 
tacts  will  exceed  "2"  so  that  chain  branching  occurs  which 
gives  rise  to  network  formation.  This  second  structural 

transition  occurs  at  the  saturation  volume  fraction  v-c  =  a' 

-  ^  c 

associated  with  the  onset  of  network  formation  (i.e.,  that 
concentration  at  which  the  average  number  of  contacts  begins 
to  exceed  "2").  From  Gurland's  observations,  the  saturation 
concentration  occurs  in  the  vicinity  of  Vf  =  0.4. 

Upon  the  onset  of  incipient  network  formation  at 
''^f  multitude  of  continuous  conduction  paths 

causes  a  transition  to  behavior  dominated  by  the  highly 
conducting  (low  resistivity)  filler  material.  The  structure 
associated  with  the  conducting  plateau  is  illustrated  in 
Figure  1  by  the  closed  circles  which  represent  contacting 
particles  participating  in  a  continuous  network  of  conduct¬ 
ing  paths.  A  further  increase  in  the  concentration  of 
conducting  particles  beyond  the  saturation  concentration 
simply  fills  the  interstices  of  the  network  with  additional 
conducting  material  so  that  the  conductivity  follows  the 
behavior  predicted  by  Equation  (3a). 

The  significance  of  network  formation  in 
establishing  high  levels  of  conductivity  was  demonstrated  by 
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Malliaris  and  Turner  [1971].  These  workers  tumbled  large 
particles  of  polymers  (polyvinyl  chloride;  polyethylene) 
with  metal  powders  (nickel;  copper).  As  illustrated  in 
Figure  2a,  the  metal  powder  adhered  to  the  surface  of  the 
larger  polymer  particles.  The  metal-powder  coated  polymer 
particles  were  subsequently  compacted  under  conditions  which 
limited  flow.  Photomicrographs  of  the  compacted  material 
revealed  metal  particle  networks  shown  in  an  idealized  form 
in  Figure  2b.  As  a  consequence  of  the  segregation  of  metal 
particles,  network  formation  (and  consequently,  high  levels 
of  conduction)  were  achieved  at  a  significantly  lower 
saturation  concentration  ((J)q  =  0.06)  than  would  be  required 
if  the  particles  v;ere  dispersed  uniformly  throughout  the 
polymer  phase. 

These  qualitative  considerations  indentify  features 
of  the  microstructure  that  must  be  incorporated  into  quanti¬ 
tative  models.  The  development  of  such  models  can  be 
approached  in  two  stages : 

o  The  development  of  transport  models  which 
account  for  a  distribution  of  path  lengths. 

o  The  development  of  statistical  models  to 
predict  path  length  distributions  in  terms 
of  the  concentration  of  filler,  and  the 
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orientation  and  shape  of  the  conducting  filler 
(e.g.,  particles,  fibers,  or  platelets). 


Approach 


Transport  Models; 

Improved  bounds  on  transport  properties  were 
recently  reported  by  Nomura  and  Chou  [1980].  These  workers 
focused  attention  on  a  system  of  aligned  fibers  of  equiva¬ 
lent  but  arbitrary  aspect  ratio,  "a".  The  results  from  this 
analysis  take  into  account  two-  and  three-point  correlation 
functions  through  the  use  of  the  modified  Eshelby  tensor 
developed  by  Wu  and  McCullough  [1977].  These  bounds  provide 
a  basis  for  establishing  the  theoretical  limits  for  conduc¬ 
tivity.  However,  without  further  treatment  these  limits  are 
of  little  practical  value  since  the  upper  and  lower  bounds 
on  electrical  conductivity  may  differ  by  14  to  16  orders  of 
magnitude.  Wu  and  McCullough  [1977]  have  shown  that 
behavior  intermediate  between  these  wide  limits  is  strongly 
dependent  upon  the  microstructure  developed  within  the 
composite  material. 


The  bounding  results  were  extended  to  account  for 
distributions  in  fiber  lengths.  Following  the  approach  of 
McGee  and  McCullough  [1981],  these  improved  bounds  were  used 
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to  obtain  model  relationships  to  predict  intermediate 
behavior  between  the  limiting  bounds  as  dictated  by  distri¬ 
butions  of  contact  path  lengths  for  collections  of 
particles,  fibers,  and  platelets.  This  model  provides  con¬ 
stitutive  relationships  for  a  "Representative  Volume 
Element"  (RVE).  The  predicted  properties  of  the  RVE  provide 
subsequent  input  to  an  "Aggregate  Model"  [Jarzebshi,  McGee 
and  McCullough,  1982],  which  accounts  for  orientation 
distributions . 

Statistical  Models; 

For  a  given  number,  N,  of  filler  elements,  the  path 
length  distributions  may  be  obtained  by  the  application  of 
"random  flight  statistics"  traditionally  used  to  predict 
dimensions  of  polymer  molecules.  The  key  parameter  of  the 
resulting  distribution  function  is  the  "contact  probabil¬ 
ity".  Studies  of  contact  probabilities  for  particles  are 
available  from  the  chemical  engineering  literature  related 
to  packed  bed  catalytic  reactors;  e.g.,  Suzuki  et  al  [1981]. 
The  results  from  these  studies  v/ere  used  to  estimate  the 
concentration  dependence  of  the  contact  probability  and 
hence  the  behavior  of  the  path  length  distributions  as  a 
function  of  particle  concentration.  These  results  were 
extended  to  fibers  through  approximate  geometrical 
constructions . 


SUMMARY  OF  MODEL  RELATIONSHIPS 


The  mathematical  structure  of  the  network  model  is 
summarized  in  Tables  I  through  III.  The  K's  represent  any 
one  of  the  following  transport  properties: 

o  electrical  conductivity 
o  thermal  conductivity 
o  diffusivity 

o  dielectric  constants  (both  real  and  imaginary) 
o  magnetic  permeability 

The  distinction  between  the  various  transport  properties  is 
achieved  by  specifying  certain  parameters  (described  in  a 
subsequent  section)  related  to  surface  and  volume  transport 
mechanisms;  however,  attention  will  be  restricted  to 
electrical  conductivity.  The  network  model  accounts  for  the 
following  important  structural  features  of  composite 
materials : 

o  intermediate  states  of  fiber  orientation 
through  an  orientation  parameter,  f  ; 
o  filler  shape  through  the  aspect  ratio 
parameter,  "a".  (For  continuous  fibers. 


20 


21 


TABLE  I 

Orientation  Dependence 

=  Bulk  transport  property 

* 

=  RVE  transport  property 
f  =  Orientation  parameter 


Axial  Orientation:  f; 


<K^>  =  i  K*  +  ^  (K*  -  Kp 


<K2>  =  <K3>  =  1  k*  +  |k*  -  i  (K*  -  K*) 


fg  =  j  C3<cos^<j)>^  -  1] 

2  2 

<cos  ({)>  =  /  cos  (|)  n  ((}))  sin4)  d^i 

a  0  P 


Planar  Orientation:  f 


<Ki>  =  1  K*  +  iK^  +  if  (K*  -  Kp 


1  *  1  * 
<Ki>  =  2  ^1  ^  2^S 


"  H 


<K3> 


=  K, 


f  =  2<cos  (1)>  -  1 

p  ^  p 

2  2 

<cos  =  /  cost})  n  (cj))  sin(j)  d(}) 

P  0  P 
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TABLE  II 

Constitutive  Properties  of  RVE 


Transport  Property 
Longitudinal 
Transverse 
Volume  Fraction 


Filler 


Kpi 

KF2 

Vf 


Matrix 


K 


m 


K 


'm 


(i  =  1,2) 

Kpi  >  Km 


K  K  {K  +  (v  -  v^)  AK  <y  >} 

*  F.rni  m  i.  li 

K  =  - i - - - - - - 

i  K_  K  +  [v^  K_  -  v^  K  ]  AK.  <y.> 
F.  m  mF.  fm 

1  1 


K  =  v^K^  +  V  K 
f  mm 


AK.  =  K„  -  K 
1  Fjj^  m 


=  <yj>  =  _  [1  -  H(Vj,  a)] 


y  (a)  =  1  -  h(a) 
o 


( continued ) 
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(Table  II,  continued) 


a 

a  >  1 


a  =  1 


0<a<l 


X 
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TABLE  III 

Path  Distributions 

Vf  =  filler  volume  fraction  a  =  filler  aspect  ratio 


Critical  Concentrations 


’  /  X -1/3 

<j)^  =  a(a) 


=  6(a) 


-1/3 


Transport  Mechanism 


Bullc;  a  =  2/3,  g  =  1/3  Surface:  a  =  0.4,  g  =  0.3 


Distribution  Parameters  (i  =  1,2) 


Mean : 


Z  .  d)  +  Z  .(b 

1  ^  C  3.  C 


Z  .  +  Z  . 
1  1 


Spread : 


6  -  (j) 

^  c  ^  c 

I  M 

Z.  +  Z. 
1  1 


(Z,)2 


(Z,)2 


=  2  In 


In 


(1-<|,  .)^R.yo 


L  3{(t)^(l-yQ)  +  (l-(t)^yQ>  J 

l-h(a)  (Table  II) 


R.  = 


/K 

F.  m 
1 


(Table  II) 


( continued ) 
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(Table  III,  continued) 


Distribution  Functions 


(i  =  1,2) 


H^(v  ,a)  =  1  -  '5(X^) 


X. 

1 


V  _  <  IJ  . 

r  1 

Hi(Vf,a)  =  $(|  xi]  ) 


<5(X^) 


e 


-T 


dT 
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a  ->■  »  7  for  short  fibers,  a  >  1  ;  for 
particles,  a  =  1  ;  and  for  platelets, 

0  <  a  <  1  .) 

These  important  features  are  discussed  in  the  following 
sections . 

Orientation  Dependence 

Under  the  Aggregate  Model  [Jarzebski,  McGee  and 
McCullough,  1982],  the  properties  of  a  bulk  specimen  may  be 
obtained  by  averaging  the  properties  of  a  transversely 
isotropic  Representative  Volume  Element  (RVE),  characterized 
by  the  properties  K*  and  K*  ,  over  appropriate  orientation 
distributions.  The  details  of  the  averaging  procedures  are 
given  elsewhere  [McGee,  1982].  The  results  for  contracted 
second  rank  tensors  corresponding  to  transport  properties, 

K,  are  summarized  in  Table  I  for  two  important  types  of 
orientation  distributions  —  axial  and  planar.  The  quanti¬ 
ties  <Ki>,  <K2>,  and  <K3>  are  the  observed  properties  of  the 
bulk  specimen  along  the  specimen  axes  (e.g.,  the  longitudi¬ 
nal  "1"  direction,  the  transverse  or  width  direction,  "2", 
and  the  perpendicular  or  thickness  direction,  "3").  The 
quantities  K*  and  are  the  properties  of  the  transversely 
isotropic  RVE  predicted  by  the  constitutive  relationship 
given  in  Table  II. 
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Axial  orientation  distributions  are  associated  with 
processing  procedures  (e.g.,  injection  molding)  in  which 
fibers  achieve  a  three-dimensional  distribution-  The  impor¬ 
tant  features  of  an  axial  distribution  are  ilustrated  in 
Figure  3.  In  this  case,  the  orientation  distribution  func¬ 
tion,  n(4)),  gives  the  fraction  of  filler  elements  which  lie 
on  a  cone  making  an  angle  (j)  with  respect  to  the  "1"  axis  of 
the  specimen.  The  normalized  distribution  is  taken  as  sym¬ 
metric  around  the  "1"  axis,  n((f))  =  n(-<}i),  and  across  the  2-3 
plane,  n((j))  =  n(iT+(|)).  The  statistical  characteristics  of 
the  axial  distribution  may  be  described  by  a  parameter,  f^, 
scaled  to  the  range  -1/2  <  fa  <  1  .  This  parameter  is 

related  to  the  second  moment  of  the  three-dimensional 
distribution : 


fa  =  ( 1/2)  [3<cos2(})>  -  1] 


(Table  I) 


with 


■it/2 


<  cos- 


■9>  = 


n  ( 4) )  cos2(j)sin(J)d.^) 


A  value  of  f^  =  -1/2  corresponds  to  a  situation 
the  fibers  are  confined  to  lie  in  the  2-3  plane 
<cos2(|)>  =  0);  a  value  of  fa  =  1  denotes  perfect 
along  the  "1"  axis  (i.e.,  <cos2,^>  =  i).  States 


in  which  all 
( i • e. , 

collimation 
of  partial 
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three-dimensional  orientation  are  indicated  by  intermediate 
values  of  ^  three-dimensional  random  orientation  is 

characterized  by  =  0  (i.e.,  <cos2<|)>  =  1/3).  In  this 
event,  the  bulk  specimen  is  isotropic  since 

<K]!>  =  <K2>  =  <K3>  .  Alternately,  an  isotropic  bulk  speci¬ 
men  is  obtained  if  the  RVE  is  isotrooic  so  that  K*  =  K*  . 

"  1  2 

In  this  situation,  the  terra  [K*  -  K*2  of  Table  I  vanishes  so 
that  the  specimen  properties  are  independent  of  f^  • 

Sheet  molding  materials  as  well  as  other  systems 
produced  by  compression  molding  can  be  characterized  by  a 
two-dimensional  planar  orientation  in  which  the  fibers  are 
confined  to  lie  in  a  plane.  In  some  cases,  injection  and/or 
transfer  molding  procedures  will  give  rise  to  stratified 
layers  that  may  be  locally  characterized  as  planar  orienta¬ 
tions.  The  features  of  this  important  orientation  distribu¬ 
tion  are  illustrated  in  Figure  4.  For  this  case,  the 
normalized  orientation  distribution,  n((|)),  gives  the 
fraction  of  fibers,  within  the  1-2  plane,  which  makes  an 
angle  with  respect  to  the  "1"  axis  of  the  specimen.  The 
distribution  is  taken  as  symmetric  across  the  1-3  plane, 
n((t))  =  n(-(t)),  and  across  the  2-3  plane,  n(.j>)  =  nC-ii+ir)  .  The 
statistical  characteristics  of  the  planar  distribution  may 
be  described  by  a  parameter,  fp  ,  scaled  to  the  range 

<  f  <  1  .  As  before,  f  is  related  to  the  second  moment 
iy  p 

of  the  two-dimensional  distribution: 


Schematic  definition  of  the  features  of  a  two-dimensional  planar 
orientation  distribution. 
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with 


fp  =  2<cos2(j)>  -  1 


(Table  I) 


f  1 

<cos2(f,>  =  J  n ( 4) )  cos2(j)d.|i  . 

0 

A  value  of  fp  =  1  corresponds  to  fiber  alignment  along  the 
"1"  axis  of  the  specimen;  a  value  of  fp  =  -i  corresponds  to 
alignment  along  the  "2"  axis.  Partial  states  of  planar 
orientation  are  described  by  intermediate  values  of  f  . 

ir 

A  value  of  fp  =  o  corresponds  to  a  quasi-isotropic  material 
with  <K]_>  =  <K2>  but  distinct  from  <K3>  .  In  the  event  that 
the  RVE  is  isotropic  such  that  K*  =  K*  ,  the  bulk  specimen 
exhibits  isotropic  character  with  <Ki>  =  <K2>  =  <K3>  . 

Recently,  McGee  [1982]  generalized  the  treatment  of 
orientation  distributions  to  account  for  distributions 
intermediate  between  three-dimensional  axial  and  two- 
dimensional  planar.  In  this  extension,  the  projection  of  a 
fiber  onto  the  1-2  plane  is  specified  by  an  angle  (|)  with  an 
associated  orientation  parameter,  fp  .  The  orientation  with 
respect  to  the  "3"  axis  is  specified  by  an  angle  6  with  an 
associated  orientation  parameter  f^  .  Under  this  revision, 
the  parameter  fg^  is  a  measure  of  the  extent  to  which  fibers 
are  tilted  out  of  the  1-2  plane  toward  the  "3"  axis.  These 
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generalized  orientation  descriptors  are  defined 
schematically  in  Figure  5. 

The  three-dimensional  axial  and  two-dimensional 
planar  distributions  are  contained  as  special  cases.  The 
two-dimensional  planar  distribution  corresponds  to 
fa  =  -1/2  .  The  three-dimensional  axial  distribution  is 
obtained  with  fp  =  o  ;  however,  in  this  revision,  the  prin¬ 
ciple  axis  of  axial  alignment  corresponds  to  the  "3"  axis. 

A  three-dimensional  random  orientation  is  obtained  with  both 
fa  and  fp  equal  to  zero. 

The  properties  of  the  bulk  specimen  are  given  in 
terms  of  the  generalized  distribution  by 

<K^>  =  +  kCfp  -  ^(1  +  fp)]  (4a) 

<K,>  =.  K  +  kCt  +  f  (1  -  f„)]  (4b) 

2  o  Pa.  P 

<K, >  =  K  +  2kf  (4c) 

3  o  a 

with 

=  (1/3)CkJ  +  2K2] 

k  =  (1/3)CK*  -  K*] 

where  K*  and  K*  are  again  the  properties  of  the  RVE  as 
predicted  by  the  constitutive  relationships  of  Table  II. 


2  <  cos^<^  >-l 
[8  <cos'^(p)-Z] 


■  '2  C0S^8  >-!] 

Qa  =  “4  [5<C0S^8>-|] 


with 

TT/Z  - 

<cos^<^>=j  r\{cp)  cos^  (p6(p 

0 

IT 

<cos"8>=j  n(8) sin8cos"8d8 

0 


Schematic  definitj.on  of  the  orientation 
parameters  of  the  generalized  orientation 
distribution.  After  McGee  ri9S2]. 
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These  important  descriptors  of  fiber  orientation 
must  be  determined  by  independent  experimental  characteriza¬ 
tions  such  as  the  laser  scattering  methods  described  by 
McGee  and  McCullough  [1982]. 

Constitutive  Relationships 

The  quantities  K*  (i  =  1,2)  of  Table  I  correspond  to 
the  properties  of  a  transversely  isotropic  Representative 
Volume  Element  (RVE)  of  the  material  system.  The  relatidn- 
ship  of  the  K*  to  the  properties,  concentration,  and  phase 
geometry  of  the  components  are  summarized  in  Table  II. 

The  quantity  Kj^  refers  to  the  appropriate  transport 
property  of  an  isotropic  matrix.  The  refer  to  the 

effective  transport  property  of  a  transversely  isotropic 
chain  of  contacting  filler  elements  with  corresponding 

to  the  longitudinal  property  and  Kp2  corresponding  to  the 
transverse  property.  In  the  event  that  the  filler  is 
isotropic,  then  Kpj^  =  Kf2  • 

The  quantities  Vf  and  vm  represent  the  volume 
fraction  of  filler  and  matrix,  respectively,  with 

Vf  -l"  V]^  —  1 . 

The  term  <yi>  is  a  structure  factor  which  takes  into 
account  the  shape  of  the  filler,  path  length  distributions, 
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and  the  underlying  mechanism  of  transport.  The  role  of 
these  various  quantities  is  further  specified  in  Table  III. 

The  quantity  y^  is  determined  by  the  true  aspect  ratio,  "a", 
of  the  filler.  The  functional  dependence  of  y^  on  the 
aspect  ratio  is  given  for  three  important  shapes  in 
Table  II;  viz.,  fibers  (a  >  1),  spherical  particles  (a  =  1), 
and  platelets  (0  <  a  <  1)  . 


The  results  summarized  in  Table  II  can  be  recast 
into  a  form  similar  to  the  Halpin-Tsai  Relationship  [Aston 
et  al,  1969];  viz.. 


with 


and 


where 


* 

K. 

1 

K 

ra 


1  + 


1  - 

XiVf 

‘^Fi 

-  K 

m 

'Si 

+  5  -K 

1  m 

1  -  • 

<1^  j_(v 

>1 

V 

1 _ t 

Fl 

v^K  +  V  K„. 

f  m  m  F 1 


(5a) 


(  5b) 


As  the  term  <yi>  approaches  the  minimum  value  of  0,  the 
reinforcing  factor  approaches  =»  so  that  the  constitutive 
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relationship  reduces  to  the  simple  form  of  Equation  (3a); 
viz .  , 


K.  =  v^K^.  +  V  K 
1  f  Fi  mm 


for  <yj^>  0 


(3a) 


Under  these  limiting  conditions,  the  highly  conducting 
filler,  (Kpj_  >>  Kjj,)  dominates  the  behavior  of  the  composite 
material . 

Alternately,  if  the  term  <yi>  approaches  unity  so 
that  the  reinforcing  factor  is  given  by  5^=0,  the  consti¬ 
tutive  relationship  takes  on  the  form  of  Equation  (3b): 


m 


K. 


K 


Fi 


K 


m 


or 


P,-  = 


v^  p  .  +  v^p^ 
r  Fi  m  m 


for  <y  =  1 


(3b) 


In  this  limit,  the  behavior  is  dominated  by  the  insulating 

matrix  (p  >>  p  )  . 

m  r  1 

For  the  special  case  of  isolated  (H]^  =  H2  =  1) 
spherical  particles  (y^  =  1/3),  the  reinforcing  factor  is 
given  by 


5l  =  52  =  2\p(vf) 


(6) 
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The  Maxwell  relationship  [Maxwell,  1892],  the  Rayleigh 
equation  [Rayleigh,  1892],  the  Kerner  equation  [Kerner, 
1956],  and  the  composites  spheres  model  [Hashin,  1970],  are 
all  equivalent  and  correspond  toii;  =  l  (i.e.,  ?  =  2)  .  The 

factor  'l^i(vf)  in  Equation  (6)  is  a  concentration-property 
dependent  correction  introduced  by  the  inclusion  of  higher 
order  correlations. 

Path  Distribution  Statistics 

The  distribution  functions,  ,  account  for  the 

transport  mechanism  and  the  path  length  distributions  in 
terras  of  the  filler  shape,  orientation,  and  concentration. 
Tentative  descriptions  of  these  functions  are  proposed  in 
Table  III.  These  functions  are  dependent  upon  the  "random 
flight  statistics"  associated  with  the  distribution  of  chain 
lengths  for  a  given  collection  of  conducting  filler 
elements.  The  key  statistical  parameter  is  the  contact 
probability  as  determined  by  the  shape,  orientation,  and 
concentration  of  the  filler  elements.  An  estimate  for  the 
contact  probability  was  obtained  from  studies  of  packed  bed 
reactors  as  reported  in  the  chemical  engineering  literature; 
e.g.,  Suzuki  et  al  [1981].  The  merged  estimates  are 
summarized  in  Table  III. 
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The  parameters  a  and  8  in  Table  III  distinguish 
between  surface  transport  (e.g.,  electrical  conductivity) 
and  volume  transport  (e.g.,  thermal  conductivity). 
Accordingly,  the  proposed  model  can  be  used  for  the  entire 
class  of  transport  properties  by  the  appropriate  distinction 
between  surface  and  volume  transport  as  reflected  by  these 
parameters.  Tentative  estimates  for  these  parameters  are 
given  in  Table  III. 

It  is  evident  that  the  application  of  the  model 
relationships  summarized  in  Tables  I,  II,  and  III  requires 
tedious  computations.  An  interactive  computer  program  has 
been  developed  to  facilitate  the  use  of  these  model  rela¬ 
tionships  as  a  screening  tool  for  ranking  candidate  con¬ 
ducting  composites.  A  user's  guide  to  the  program,  with 
illustrative  examples,  is  given  in  the  Appendix. 

Interactive  computer  programs  are  available  for  predicting 
the  associated  thermoelastic  properties.  The  model  develop¬ 
ment  for  the  thermoelastic  properties  is  described  elsewhere 
[Jarzebski,  McGee,  and  McCullough,  1982]. 

Models  for  Three  Component  Composites 

As  shown  by  McCullough  and  Peterson  [1977], 
combinations  of  two  or  more  filler  materials  in  a  composite 
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system  introduce  additional  degrees  of  freedom  for  tailoring 
composite  materials  to  multiple  design  requirements. 
Accordingly,  it  is  useful  to  extend  the  preceding  model  form¬ 
ulation  to  include  systems  comprised  of  two  distinct  filler 
materials  that  may  be  differentiated  by  their  (i)  proper¬ 
ties,  (ii)  shape,  (iii)  orientation,  and/or  any  combination 
of  these  features. 

Unfortunately,  rigorous  theoretical  treatments  which 
take  into  account  higher  order  correlations  between 
different  types  of  fillers  are  not  currently  available. 
Consequently,  approximate  form.ulations  are  required  for  an 
extension  to  three  component  material  systems. 

As  a  first  order  approximation,  the  second  filler 
phase  will  be  viewed  as  an  excluded  volume.  The  exclusion 
of  a  filler  of  type  "1"  from  the  volume  occupied  by  the 
filler  phase  of  type  "2"  forces  the  conducting  elements  of 
the  type  "1"  into  a  smaller  relative  volume  and  hence  enhan¬ 
ces  chain  formation,  as  illustrated  in  Figure  2.  The 
apparent  volume  fraction  of  type  "1"  filler  in  the  reduced 
volume  of  an  intermediate  two-component  system  (System  I)  is 
given  by 


V 


I 

FI 


^Fi/d 


V 


F2 


) 


(7a) 


40 


where  vp]_  is  the  true  volume  fraction  of  the  type  "1"  filler 
and  Vp2  is  the  true  volume  fraction  of  the  filler  viewed  as 
the  excluded  volume.  The  apparent  volume  fraction  of  the 
matrix  in  System  I  is  taken  as 


{7b) 


The  transport  properties  <K^>  (i  =  1,2,3)  of  the  inter¬ 
mediate  two-component  system  may  be  computed  by  the  rela¬ 
tionships  given  in  Tables  I  -  III  and  implemented  in  the 
computational  algorithm  described  in  the  Appendix.  In  this 
application,  the  concentration  is  taken  as  the  apparent 
volume  fractions,  v^^  and  v^  ,  given  by  Equations  (7).  The 
filler  properties,  aspect  ratio,  and  orientation  parameters 
are  specified  as  those  associated  with  the  filler  of  type 
"  1  . 

A  similar  excluded  volume  argument  may  be  applied 
for  fillers  of  type  "2"  segregated  by  exclusion  from  the 
volume  occupied  by  fillers  of  type  "1".  The  apparent  volume 
fraction  of  type  "2"  filler  in  the  reduced  volume  of  the 
second  two-component  system  (System  II)  is  given  by 


with  the  matrix  volume  fraction  in  System  II  given  by 
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As  before,  the  quantities  vpi  and  vf2  are  the  true  volume 
fractions  of  filler  "1"  and  "2",  respectively. 


The  transport  properties,  <KlI>  (i  =  1,2,3),  of 
the  second  intermediate  two-component  system  may  be  computed 
by  the  relationships  given  in  Tables  I  -  III  and  implemented 
in  the  computer  programs  described  in  the  Appendix.  The 
concentrations  of  the  filler  and  matrix  components  are  taken 

as  the  apparent  volum.e  fractions,  vll  and  vH  determined 

1 2  m 

from  Equations  (3);  the  filler  properties,  aspect  ratio,  and 
orientation  parameters  are  specified  as  those  associated 
with  the  filler  of  type  "2". 


The  transport  properties,  <Kj_>  (i  =  1,2,3),  of  the 
three  component  system  may  be  estimated  as  a  combination  of 
the  properties  of  System  I  and  System  II  weighted  by  their 
respective  volume  fraction  contributions  to  the  combined 
system.,  viz.,  v^  =  1  _  v-p2  and  vjj  =  1  -  vp]_  .  However, 
the  volume  and  the  properties  of  the  matrix  phase  were  uti¬ 
lized  in  both  intermediate  systems.  Consequently,  the 
combined  system  must  be  corrected  to  account  for  this  dual 
contribution  of  the  matrix  material r  viz.. 


<K  >  =  (1  -  V  )  <K>  +  (1  -  v^^  )  <K,  >  -  V  K 

1  F2i  i'll  mm 


where  i  =  1,2,3  for  the  properties  along  the  1,  2,  and  3 
axes  of  the  specimen. 

In  the  event  that  both  filler  types  are  particulate 
(a  =  1)  and  the  concentrations  are  well  above  the  critical 
concentration  for  network  formation.  Equation  (9)  simplifies 
to 


<K> 


-  ^F1  ^F1 


^F2"^F2 


V  K 
m  m 


(10) 


The  transport  properties  computed  via  Equation  (9) 
provide  conservative  estimates  for  thermal  and  electrical 
conductivities.  The  excluded  volume  argument  ignores 
contracts  (and  higher  order  correlations)  between  the  two 
types  of  fillers.  This  limitation  may  not  be  particularly 
serious  if  the  filler  types  differ  by  several  orders  of 
magnitude  in  their  conductivities  (i.e.,  one  filler  appears 
to  be  an  insulator  relative  to  the  other) .  A  system 
comprised  of  glass  (or  graphite)  fibers  and  metal  fillers 
(particulate,  fiber,  or  platelet)  should  approach  this 


condition . 


SPECIFICATION  OF  MODEL  PARAMETERS 


The  implementation  of  the  relationships  summarized 
in  Tables  I  through  III  requires  the  following  input 
information ; 

o  the  effective  properties  of  the  filler  (Kp) 
and  the  polymer  matrix  (Kn,); 

o  the  volume  fraction  of  the  filler  (vp)  and 
resin 

o  the  state  of  orientation  as  characterized  by 
fp  and/or  fa; 

o  the  filler  geometry  as  characterized  by  the 
aspect  ratio  (a);  with  a  >  1  for  fibers, 
a  =  1  for  spherical  (or  near  spherical) 
particles,  and  a  <  1  for  platelets. 

Methods  for  assessing  these  quantities  are  described  in 
the  remainder  of  this  section. 

Effective  Properties 

The  electrical  resistivity  (p  =  l/K)  and  density  of 
several  candidate  filler  and  resin  materials  are  summarized 
in  Table  IV.  These  values  correspond  to  typical  handbook 
listings  for  the  various  pure  materials.  The  conductivity 
of  the  pure  polymeric  material  serves  as  the  appropriate 
input  parameter  (Kii;i)  for  the  continuous  binder  phase. 
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TABLE  IV 

Typical  Properties  of  the  Components 
of  Conducting  Composites 


Material 

Resistivity 
( ohm-cra) 

Density 
( gm./ cm3 ) 

A1 uminum 

2.63 

Tr 

O 

2.77 

Copper 

1.72 

X 

10"^ 

8.96 

Graphite  Fiber 

1.50 

X 

10“^ 

1 .77 

Iron 

10.0 

X 

10~® 

7.15 

Nickel 

7.8 

X 

10"^ 

8.33 

Silver 

1.47 

X 

10“® 

10.6 

Steel 

81.0 

X 

10-6 

7.78 

Bakelite  (phenolic) 

4  X 

10 

11 

1.43 

Nylon  6-10 

5  X 

10 

14 

1.08 

Polyester 

4  X 

10 

16 

1.4 

Polyphenylene  Sulfide 

1  X 

10 

16 

1.34 
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As  noted  by  Malliaris  and  Turner  [1971],  as  well  as 
Li  et  al  [1982],  the  effective  conductivity  of  chains  of 
contacting  filler,  Kp  ,  may  be  substantially  reduced  from 

that  of  the  pure  material,  K®  ;  viz., 

F 


where  is  the  conductivity  of  the  pure  filler  material  and 
and  X  is  a  reduction  factor  attributed  to  geometrical  con¬ 
tact  irregularities  associated  with  surface  roughness  and/or 
surface  coatings .( e.g. ,  oxide  layers  or  special  sizing 
agents).  Malliaris  and  Turner  [1971]  report  a  reduction 
factor  of  X  =  10  for  contacting  chains  of  nickel 
particles;  Gurland's  results  for  contacting  silver  particles 
are  consistent  with  this  reduction  factor.  Li,  Streider  and 
Joy  [1982]  report  a  reduction  factor  of  x  =  3  x  10  for 
contacting  graphite  fibers.  These  results  imply  that  the 

reduction  factor,  x  >  niay  be  taken  as  essentially  constant 

—  3 

and  on  the  order  of  1  -  3  x  10  .  The  conductivity  reduc¬ 

tion  factor  corresponds  to  an  enhancement  of  the  effective 
resistivity  by  a  factor  of  3  -  10  x  10^  . 

Concentration  of  Components 

As  a  consequence  of  volume  averaging,  volume 
fractions  emerge  as  the  natural  descriptors  of  concentration. 
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However,  weight  fraction  concentrations  are  commonly  employed 
to  specify  the  composition  of  composite  materials-  These 
measures  of  concentration  are  related  through  the  densities 
of  the  components;  viz.. 


Wf  + 


(1  -  W^) 


(12a) 


where  v/f  is  the  weight  fraction  of  filler  and  "r"  is  the 
ratio  of  the  density  of  filler  to  the  density  of  resin; 


r  =  df/dm  (12b) 

The  term  vf  is  the  volume  fraction  of  filler  with 
Vm  =  1  -  vf  . 

The  relationship  of  volume  fraction  to  weight 
fraction  of  filler  for  several  important  material  systems  is 
displayed  in  Figure  6. 


Metal  coated  fibers  are  of  particular  interest  as 
candidate  fillers  for  conducting  composites.  The  density  of 
metal  coated  fibers  is  dependent  upon  the  thickness  of  the 
metal  layer.  Accordingly,  the  relationship  between  volume 
fraction  and  v/eight  fraction  is  dependent  upon  the  relative 
thickness  of  the  coating  on  the  fiber;  viz.. 


V  „  = 


w  ^d 
f  m 


w-d  + 
f  m 


1  -  w^)  td,^  +  (d^  -  d,^^)x(2  -  A)!! 


(13a) 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


Wf 


Relationship  of  volume  fraction  filler  concentra¬ 
tion,  Vf  ^  to  weight  fraction,  wf  ,  for  several 
material  systems:  polyphenylene  sulfide/graphite 
(r  =  1.32);  Nylon  6-lO/graphite  (r  =  1.64);  poly¬ 
phenylene  sulf ide/steel  (r  =  5.31);  Nylon  6-10/ 
steel  ( r  =  7.20). 
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where  dj^  is  the  density  of  the  polymer  matrix,  dh  is  the 
density  of  the  base  fiber,  and  d^-,  is  the  density  of  the 
metal  coating.  The  parameter  "X"  is  the  ratio  of  the 
thickness  of  the  coating  ,  t^  ,  to  the  total  radius,  R  ,  of 
the  coated  fiber; 

^  =  tc/R  (13b) 

The  relationship  of  volume  fraction  to  weight  fraction  for 
graphite  fibers  coated  with  various  thicknesses  of  nickel 
are  displayed  in  Figure  7  for  a  Nylon  6-10  matrix  and  in 
Figure  8  for  a  polyphenylene  sulfide  matrix. 

State  of  Orientation 

Information  concerning  the  orientation  distribution 
of  fibers  can  be  obtained  from  photomicrographs  of  certain 
sections  of  the  test  specimen.  Sections  taken  from  a  side 
surface  and  an  end  surface  are  required  to  characterize  the 
out-of-plane  tilting  of  fibers  as  measured  by  fg^  _  Sections 

taken  from  a  top  surface  and  a  mid-plane  surface  are 
required  to  characterize  the  uniformity  of  the  in-plane 
orientation  measured  by  fp  .  These  sections  are 
schematically  defined  in  Figure  9. 

Typical  photomicrographs  of  these  four  sections  are 
shown  in  Figure  10  for  a  test  specimen  of  graphite  fibers 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

Wf 


Figure  7.  Relationship  of  volume  fraction  filler  concentra 
tion,  Vf  ,  to  weight  fraction,  wf  ,  for  graphi 
fibers  coated  with  various  thicknesses  of  nick 
in  a  Nylon  6-10  matrix.  The  parameter  \  denotes 
the  ratio  of  the  thickness  of  metal  to  the  total 
radius  of  the  fiber. 


4->  0) 


Relationship  of  volume  fraction  filler  concentra^ 
tion,  V£  ,  to  weight  fraction,  Wf  ,  for 
fibers  coated  with  various  thicknesses  o 
in  a  polyphenylene  sulfide  matrix.  The  parameter 
X  denotes  the  ratio  of  the  thickness  of  metal  to 
the  total  radius  of  the  fiber. 


Figure  8 . 
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Figure  9. 


Schematic  diagram  illustrating  various  locations 
from  which  specimen  sections  were  sampled. 
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(20%  by  weight)  in  a  polyphenylene  sulfide  matrix.  Photo¬ 
micrographs  of  four  sections  taken  from  a  test  specimen  of 
nickel  coated  graphite  fibers  (30%  by  V'/eight)  in  Nylon  6-10 
are  displayed  in  Figure  11. 

From  these  direct  observations,  the  number  of 
fibers,  N(ij;)  which  are  oriented  between  \p  ~  Af  and 
ip  +  Ail;(iJ;  =  ())  or  0)  with  respect  to  the  appropriate  specimen 
axis  could  be  counted  and  a  histogram  of  the  orientation 
distribution  constructed.  The  associated  orientation 
parameters,  f^  or  fp  ,  could  then  be  obtained  by  the 
numerical  integrations  prescribed  in  Table  I . 

Recently,  McGee  and  McCullough  [1982]  have  shown 
that  photographic  negatives  of  sample  sections  can  be  used 
as  diffraction  masks  in  a  Fraunhoffer  diffractometer  to 
obtain  diffraction  patterns  characteristic  of  the  state  of 
orientation.  The  components  of  a  Fraunhoffer  diffractometer 
are  shown  in  Figure  12.  Diffraction  patterns  obtained  from 
the  photomicrographs  shown  in  Figure  10  are  displayed  in 
Figure  13;  diffraction  patterns  from  the  photomicrographs 
given  in  Figure  11  are  shov/n  in  Figure  14. 

McGee  and  McCullough  [1982]  have  shown  that  the 
relative  intensity  of  the  diffraction  pattern  is  propor¬ 
tional  to  the  orientation  distribution  shifted  by  n/2;  viz.. 


c  o  c 

0  (0  "H 
Hlw  (0 
>iU  -P 
la  3  ^ 
03  o 


Figure  12.  Schematic  diagram  illustrating  the  components  of  a  Fraunhoffer 
diffractometer. 
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Figure  14.  Diffraction  patterns  obtained  from  the  sample  sections  of  Figure  11 
for  a  30%  (by  weight)  Nylon  6-lO/nickel  coated  graphite  fiber 
specimen.  After  McGee  [1982]. 
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I(r,>|;)  =  I(r)n(iJ)  +  'ir/2)  (14) 

where  r  and  i};  are  the  polar  coordinates  of  the  diffraction 
pattern  as  measured  from  the  central  spot  as  the  origin. 
Accordingly,  an  intensity  trace  over  a  circle  of  radius  "r" 
on  the  diffraction  pattern  yields  the  orientation  distribu¬ 
tion  function.  Optical  densitometer  traces  of  the  diffrac¬ 
tion  patterns  shown  in  Figures  13  and  14  are  displayed  in 
Figures  15  and  16.  Procedures  for  extracting  orientation 
parameters  from  these  intensity  profiles  are  given  elsewhere 
[McGee,  1982]  .  Application  of  these  procedures  shows  that 
both  the  graphite  fiber/polyphenylene  sulfide  and  the  nickel 
coated  graphite/Nylon  6-10  specimens  are  in  a  planar  state 
of  orientation  ( f ^  ^  -  1/2).  The  graphite  fiber/polypheny- 
lene  sulfide  specimen  shown  in  Figure  10  exhibits  a  high 
degree  of  planar  orientation  with  f  s  0.7  .  In  contrast, 
the  nickel  coated  graphite  fiber/Nylon  6-10  specimen  shown 
in  Figure  11  exhibits  a  nearly  in-plane  random  orientation 
with  fp  =  0.2-0. 3.  It  will  be  shown  in  the  subsequent 
section  that  these  different  states  of  orientation  are 
reflected  in  the  observed  resistivities  of  the  specimens. 

Fiber  Aspect  Ratio 

Estimates  of  the  average  aspect  ratio  can  be 
obtained  from  photomicrographs  such  as  those  displayed  in 
Figures  10  and  11.  Direct  measurements  of  the  ratio  of 


Optical  densitometer  tracings  of  the  diffraction  patt 
Figure  13  obtained  from  sections  of  a  20%  (by  weight) 
sulfide/graphite  fiber  specimen.  After  McGee  [1982]. 
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fiber  lengths  to  fiber  widths  yield  aspect  ratios  on  the 
order  of  20-30  for  the  graphite  fiber/polyphenylene  sulfide 
specimen  (Figure  10)  and  aspect  ratios  on  the  order  of  30-50 
for  the  nickel  coated  graphite  fiber/Nylon  6-10  specimen 
( Figure  11 )  . 

Alternately,  information  concerning  the  aspect  ratio 
can  be  obtained  from  diffraction  patterns  such  as  those 
shown  in  Figures  13  and  14.  Recently,  McGee  [1982]  has 
reported  an  analysis  of  the  radial  intensity  profiles  of 
diffraction  patterns  that  yields  the  average  aspect  ratio  of 
■(-he  fibers.  It  is  anticipated  that  future  developments  in 
the  analyses  of  such  diffraction  patterns  wixl  yield  guanti— 
^^(-ive  information  concerning  aspect  ratio  distr ioutions 
and/or  particle  size  distributions. 


COMPARISONS  OF  PREDICTIONS  WITH  EXPERIMENTAL  DATA 


The  resort  to  a  modeling  approach  requires  that  the 
suitability  of  the  simplifying  abstractions  be  tested  by 
comparisons  between  predicted  and  measured  properties.  The 
relationships  summarized  in  Tables  I  through  III  point  out 
the  necessity  for  conducting  independent  experimental 
characterizations  to  determine  the  crucial  descriptors  of 
the  microstructure;  the  orientation  parameters,  fp  and/or 
fa  ,  and  the  aspect  ratio.  Unfortunately,  such  structural 
information  is  available  for  only  a  limited  number  of 
systems  reported  in  the  literature.  Particulate  filled 
systems  provide  the  simplest  system  for  a  test  of  the  model 
relationships.  For  these  isotropic  materials,  the  bulk 
properties  are  independent  of  orientation  and  the  aspect 
ratio  may  be  taken  as  unity.  Comparison  of  measured  and 
predicted  properties  of  short-fiber  reinforced  systems 
requires  estimates  of  the  state  of  orientation  and  the 
aspect  ratio. 

Particulate  Systems 

Gurland  [1966]  reported  electrical  resistivities  of 
silver  particles  in  a  phenolic  resin  for  particulate 
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concentrations  ranging  from  10  to  50%  (by  volume) .  The 
silver  particles  v^ere  essentially  of  the  same  size.  A  com¬ 
parison  of  the  model  predictions  with  Gurland's  measurements 
is  shown  in  Figure  17. 

The  model  predictions  are  in  good  agreement  with  the 
resistivity-concentration  profiles  observed  by  Gurland.  In 
particular,  the  sharp  transition  in  resistivity  is 
accurately  predicted.  A  comparison  of  the  data  reported  by 
Kv/an  et  al  [1980]  for  silver  coated  glass  spheres  in 
polyester  shov/ed  good  correspondence  with  predicted  values. 
Improved  agreement  could  be  achieved  by  slight  modifications 
of  the  values  for  the  transport  mechanism  parameters  (a.g., 
a  =  0.35,  B  =  0.25) . 

Resistivity  data  reported  by  Aharoni  [1977]  for  iron 
particles  in  a  polyimide-amide  matrix  could  not  be  recon¬ 
ciled  with  the  model  predictions.  Aharoni ' s  results  indi¬ 
cated  a  transition  from  insulating  to  conducting  behavior 
for  volume  fractions  in  the  vicinity  of  =  0.1.  Nicodemo 
et  al  [1978]  attribute  this  anomalously  low  transition  con¬ 
centration  to  nonuniforraity  in  the  size  distribution  of  the 
iron  particles.  These  workers  propose  that  a  broad  distri¬ 
bution  of  particle  size  enhances  particle-particle  contacts. 
Such  enhanced  contact  probabilities  would  decrease  the  con- 
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Figure  17.  Comparison  of  model  predictions  with  resistivity  data  for  sil 
particles  in  a  Bakelite  (phenolic)  matrix.  Predictions  are  g 
by  the  solid  line;  resistivity  data  reported  by  Gurland  [1966 
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centration  of  particles  required  for  chain  and  network 
formation • 

Short-Fiber  Systems 

The  influence  of  fiber  aspect  ratio  on  the 
resistivity-concentration  profile  is  illustrated  in 
Figure  18  for  fibers  in  a  three-dimensional  random  orien¬ 
tation.  The  profiles  are  based  on  silver  fibers  in  a  poly¬ 
propylene  matrix.  This  result  predicts  that  the  critical 
volume  fraction  with  increasing  aspect  ratio  is  a 
consequence  of  enhanced  contact  between  longer  fibers  in  a 
random  state  of  orientation. 

The  decrease  in  the  critical  concentration  with 
increasing  aspect  ratio  illustrated  in  Figure  13  is  in  good 
agreement  with  observations  reported  by  Bradbury  and  Biggs 
[1980].  These  workers  reported  limited  resistivity- 
concentration  profiles  for  aluminum  fibers  (aspect  ratio 
12.5  and  24)  in  a  polypropylene  m.atrix.  The  state  of  orien 
tation  of  the  fibers  was  not  specified.  If  the  state  of 
fiber  orientation  is  assumed  to  be  three-dimensional  random 
good  agreement  is  achieved  between  the  predicted  and 
observed  conductivities . 

Unpublished  resistivity  measurements  for  several 
short-fiber  systems  was  provided  by  Trabacco  [1982].  The 
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molded  test  specimens  were  treated  by  three  different 
methods  to  achieve  contact  between  the  specimens  and  elec¬ 
trodes;  (i)  the  specimens  were  painted  with  a  conducting 
paint,  (ii)  the  specimens  were  plated  with  a  conducting 
metal,  and  (iii)  the  specimens  were  plated  and  dessicated. 

Selected  specimens  were  sectioned  and  subjected  to 
the  orientation  analyses  described  in  the  previous  section. 
These  analyses  revealed  that  the  state  of  fiber  orientation 
varied  from  specimen  to  specimen.  This  variation  was 
attributed  to  processing  induced  variations  in  the  micro¬ 
structure.  The  structural  variations,  in  turn,  lead  to 
significant  data  scatter  between  specimens.  These  findings 
implied  that  detailed  structural  analyses  should  be  con¬ 
ducted  on  each  specimen  in  order  to  assess  the  inherent 
conductivity  of  the  candidate  material  system.  In  order  to 
circumvent  the  necessity  for  such  a  detailed  characteriza¬ 
tion  program,  an  alternate  data  analysis  was  pursued  that 
permitted  the  utilization  of  the  existing  data. 

Inspection  of  the  photomicrographs  of  Figures  10  and 
11  shows  that  the  fiber  aspect  ratios  of  these  samples  are 
on  the  order  of  20  (or  more).  The  predictions  displayed  in 
Figure  18  suggest  that  short-fiber  systems  with  this  aspect 
ratio  are  within  the  conducting  plateau  for  volume  fraction 
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filament  geometries  are  indicated  by  a  =  1  through  a  =  20 
respectively.  The  profiles  are  based  on  silver  filler  in 
dimensional  random  state  of  orientation  in  a  polypropylen 
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concentrations  in  excess  of  0.1.  In  this  regime,  the 
behavior  is  dominated  by  the  conducting  filler  (Kp  >>  Km)  so 
that  the  specimen  conductivity  can  be  reasonably  approxi¬ 
mated  by  a  simplified  form  of  Equations  (4); 

<K^>  .  Kj,(l  +  fp)  (1  -  f^) 

<K2>  =  Kpd  -  fp)  (1  -  f^) 

<K3>  =  Kp(l  .  2f^) 

with 

Kj,  =  (l/3)v^Kp 

where  Kp  is  the  effective  conductivity  of  chains  of 
contacting  filler  and  is  related  to  the  conductivity  of  the 
pure  filler  through  Equation  (11).  This  formulation  takes 
into  account  th,e  possibility  for  out-of-plane  tilting  of  the 
fibers.  The  positive  deviations  of  the  parameter,  f^  ,  from 
-0.5  serve  as  a  measure  of  out-of-plane  fiber  tilt  and,  con¬ 
sequently,  a  measure  of  the  deviation  from  a  state  of  planar 
orientation . 

These  simplified  relationships  provide  three 
equations  with  three  unknown  parameters  —  /  fp  /  and 

fa  .  These  parameters  can  be  computed  from  the  measured 
resistivities  by  the  following  relationships: 


(15a) 

(15b) 

(15c) 

(15d) 
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Kp  =  [<K^>  +  <K2>  +  <K3>]/v^  (16a) 

fp  =  C<K^>  -  <K2>]/C<K^>  +  <K2>]  (16b) 

2<K2>  -  C<K^>  +  <K^>2 

"a  ”  2[<K^>  +  <K^>  +  <K3>]  (16c) 

with  <K^>  =  l/<p^>  where  <p^>,  <P2^ /  and  <p3>  are  the 
measured  resistivities  in  the  longitudinal,  width,  and 
thickness  directions,  respectively. 

The  parameters  computed  from  the  data  reductions 
of  Equations  (16)  provide  a  means  for  isolating  the  inherent 
characteristics  of  each  test  specimen.  These  inherent 
characteristics  provide  a  rational  basis  for  identifying  the 
origins  of  data  scatter  between  various  specimens.  This 
analysis  was  applied  to  three  candidate  material  systems: 
polyphenylene  sulfide/graphite  fiber;  Nylon  6-lO/nickel 
coated  graphite  fiber;  and  Nylon  6-10/steel  fiber. 

Polyphenylene  sulfide/graphite  fiber: 

Resistivity  data  for  polyphenylene  sulfide  filled 
with  graphite  fibers  [Trabacco,  19S2]  are  summarized  in 
Table  V  and  displayed  in  Figure  19.  The  display  of  the  data 
in  Figure  19  shows  that  the  scatter  between  measurements 
i^^de  on  similarly  treated  samples  is  essentially  equivalent 


Resistivity  Data  for  Polyphenylene  Sulfide/Graphite  Fiber  Composites 
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Figure  19.  Resistivity  data  for  polyphenylene  sulfide/ 
graphite  fiber  composites  [Trabacco,  1982]. 
Longitudinal  resistivities  are  given  by  O  ; 
transverse  (width)  resistivities  are  indicated 
;  perpendicular  (thickness)  resistivities 
are  given  by  A  .  Open  symbols  denote  specimens 
plated  and  dessicated?  closed  symbols  denote 
plated  specimens;  painted  specimens  are  identi¬ 
fied  by  the  partial  closure  of  the  symbols. 
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to  the  variation  of  the  measurements  taken  on  specimens  sub¬ 
jected  to  the  different  treatments  —  (i)  painted,  (ii) 
plated,  and  (iii)  plated  and  dessicated.  Accordingly,  the 
data  from  the  various  treatments  may  be  pooled.  The  data 
display  further  shows  the  expected  trend  for  the  longitudi¬ 
nal,  <pi>,  and  width,  <P2> »  resistivities  to  decrease  with 
increasing  concentation  of  the  conducting  graphite  filler. 

The  inherent  parameters  computed  by  Equations  (16) 
are  summarized  in  Table  VI.  This  data  reduction  shows  that 
all  specimens  were  essentially  in  a  planar  state  of  orien¬ 
tation  since  all  f^  values  were  very  near  -0.5.  This  result 
is  consistent  with  the  direct  orientation  measurements  shown 
in  Figures  10,  13,  and  15.  Most  of  the  specimens  of  the  20% 
(by  weight)  material  tend  to  achieve  a  higher  state  of 
in-plane  orientation  (f  =  0.6-0. 8)  than  the  specimens  of  . 

Ir 

the  30%  (by  weight)  material.  Again,  these  results  are  con¬ 
sistent  with  the  direct  orientation  measurements  illustrated 
in  Figures  10,  13  and  15  for  a  specimen  of  the  20%  (by 
weight)  material.  Most  of  the  specimens  of  the  30%  (by 
weight)  material  tend  toward  an  in-plane  random  state  of 

orientation  (f  =  0.2-0. 4)  so  that  the  longitudinal  and 
P 

width  resistivities  are  essentially  equivalent.  (The  nega¬ 
tive  sign  on  some  of  the  f^  values  implies  an  interchange  or 
mislabeling  of  the  "1"  and  "2"  indices  used  to  identify  the 
length  and  width  directions.) 


Computed  Parameters  for  Polyphenylene/Graphite  Fiber  Composites 
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The  effective  chain  resistivity,  relatively 

constant  with  an  average  value  of  0.47  ohm-cm.  This  value 
of  the  effective  resistivity  is  in  good  agreement  with  the 
resistivity  enhancement  factor  (1/x  =  3  x  102) ,  reported  by 
Li  et  al  [1982],  applied  to  the  resistivity  of  pure  graphite 
fibers  (Pp  =  1*5  x  10~^  ohm-cm). 

This  analysis  indicates  that  the  data  scatter 
shown  in  Figure  19  is  a  consequence  of  processing  induced 
variations  in  the  state  of  fiber  orientation  rather  than 
measuring  errors. 

Nylon  6-10/nickel  coated  graphite  fiber; 

Resistivity  data  for  Nylon  6-10  filled  with  nickel 
coated  graphite  fibers  are  summarized  in  Table  VII 
[Trabacco,  1982]  and  displayed  in  Figure  20.  Examinations 
of  photomicrographs  show  that  the  thickness  of  the  nickel 
coating  corresponds  to  10  to  20%  of  the  fiber  radius  (i.e., 

X  =  0.10-0.20);  accordingly,  the  volume  fraction  of  filler 
may  be  obtained  from  Equations  (13)  or  Figure  7..  The  data 
displayed  in  Figure  20  shows  the  remarkable  trend  for  the 
resistivity  to  increase  as  the  concentration  of  the 
conducting  filler  is  increased. 

Data  reductions  obtained  from  the  application  of 
Equations  (16)  are  summarized  in  Table  VIII.  This  analysis 


TABLE  VII 
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Obtained  from  Figure  7  v/ith  X  =  0.10  -  0.20 


Resistivity,  ohm  cm 


10 
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Iwl 


ra«i 


Resistivity  data  for  Nylon  6-10/nickel  coated 
graphite  fiber  composites  [Trabacco,  1982]. 
Longitudinal  resistivities  are  given  by  0  ; 
transverse  (width)  resistivities  are  indicated 
;  perpendicular  (thickness)  resistivities 
are  given  by  A  .  Open  symbols  denote  specimens 
plated  and  dessicated;  closed  symbols  denote 
plated  specimens;  painted  specimens  are  identi¬ 
fied  by  the  partial  closure  of  the  symbols. 


Computed  Parameters  for  Nickel  Coated  Graphite  Fiber/Nylon  6-10 
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again  shows  that  all  the  specimens  of  the  Nylon  6-lO/nickel 
coated  graphite  fiber  materials  were  in  a  planar  state  of 
orientation  since  f^  =  -0.5.  This  result  is  consistent  with 
the  direct  orientation  measurements  shown  in  Figures  11,  14 
and  16.  In  contrast  to  the  polyphenylene/graphite  fiber 
material,  the  lower  concentrations  (30%  by  weight,  10%  by 
volume)  of  the  Nylon  S-lO/nickel  coated  graphite  fiber 
materials  tended  toward  a  random  in— plane  orientation 

^  ^p  -  0*3)  while  the  higher  concentration  specimens  (40%  by 
weight;  15%  by  volume)  exhibited  a  higher  degree  of  in-plane 
alignment  ( f  ^  =  0.7-0. 8).  Again,  these  results  are  consis- 
with  the  direct  orientation  measurements  shown  in 
Figures  11,  14,  and  16  for  a  30%  (by  weight)  specimen. 

The  effective  resistivity  of  the  30%  (by  weight) 
material  was  essentially  constant  at  5  x  10-4  ohm-cm.  Con¬ 
sistent  with  the  trends  of  Figure  20,  the  average  effective 
resistivity  of  the  40%  (by  weight)  material  increased  to 
25  X  10-4  .  This  five-fold  increase  in  effective  resis¬ 
tivity  is  due  to  degradation  of  the  metal  coating.  A  typi¬ 
cal  photomicrograph  of  a  40%  (by  weight)  specimen  is  shown 
in  Figure  21  [Shaffer,  1982].  The  fibers  were  exposed  by 
disolving  the  Nylon  matrix  with  phenol.  This  photomicro¬ 
graph  shows  that  the  metal  coatings  have  undergone  dedrada- 
Evidently,  sufficient  fiber-fiber  abrasion  occurs  at 


tion . 


2000x 


Figure  21.  Scanning  electron  micrograph  of  a  40%  (by 
weight)  specimen  of  nickel  coated  graphite 
fibers  from  a  Nylon  6-10  matrix  at  2000x 
[Shaffer,  1982].  The  Nylon  6-10  matrix  was 
removed  with  phenol. 
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these  higher  volume  fraction  loadings  to  damage  the  fragile 
metal  coatings. 

This  observation  of  fiber  damage  at  the  higher 
concentrations  coupled  with  the  fact  that  a  10%  increase  in 
weight  fraction  yields  only  a  5%  increase  in  volume  fraction 
(see  Table  VII)  strongly  suggests  that  conducting  composites 
comprised  of  metal  coated  fibers  should  be  restricted  to 
relatively  low  concentrations  of  the  coated  fibers. 

• 

Nylon  6-10/steel  fiber; 

Resistivity  data  for  Nylon  6-10  filled  with 
stainless  steel  fibers  [Trabacco,  1982]  is  summarized  in 
Table  IX.  These  specimens  correspond  to  a  volume  fraction 
fiber  of  0.01.  The  data  reductions  obtained  through  the 
application  of  Equations  (16)  are  also  summarized  in 
Table  IX.  As  was  the  case  for  the  preceding  specimens,  the 
Nylon  6-lO/steel  fiber  material  was  in  a  state  of  planar 
orientation  with  f^  =  -0.49.  This  material  exhibited  a  ten¬ 
dency  to  achieve  a  moderate  state  of  in-plane  orientation 
with  f^  =  0.4-0. 6.  The  effective  resistivity  was 
0.8  X  10“3.  This  value  is  consistent  with  the  standard 
reduction  factor  (1/x  =  10-3)  applied  to  the  resistivity  of 
steel  listed  in  Table  IV. 

These  results  suggest  that  the  long  sturdy  steel 
filaments  withstood  processing  degradations  to  yield 
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respectable  conductivity  (<pi>  =  0.3)  at  very  low  volume 
fractions  (vp  =  0.01). 

Summary 


These  comparisons  of  the  model  predictions  with 
experimental  observations  on  particulate  and  short-fiber 
systems  shown  that  the  model  formulation  provides  a  rational 
basis  for  explaining  the  behavior  of  conducting  composites. 
The  agreement  between  predicted  and  measured  values  lends 
support  to  the  use  of  the  model  relations  as  a  screening 
tool  for  selecting  candidate  materials.  Furthermore,  the 
consistency  of  the  computed  orientation  parameters  with 
direct  orientation  measurements  lends  confidence  to  the  use 
of  the  modeling  relationships  to  explore  the  influence  of 
processing  induced  changes  in  the  state  of  orientation  and 
aspect  ratio  on  achievable  conductivities . 

The  analyses  of  the  various  short-fiber  systems 
emphasizes  the  important  role  that  processing  induced 
variations  in  the  microstructure  play  in  establishing  bulk 
conductivity.  The  dependence  of  bulk  conductivity  on  pro¬ 
cessing  induced  variants  of  the  microstructure  points  out 
the  futility  of  attempting  to  compare  properties  of 
specimens  of  unknown  microstructure  and  underscores  the 
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importance  of  conducting  structural  characterizations  in 
conjunction  with  traditional  testing  programs. 


CONCLUSION 


The  incorporation  of  the  concept  of  chain  and 
network  formation  of  contacting  conductors  into  a  model  for 
electrical  transport  provides  a  rational  basis  for 
explaining  sharp  transitions  from  insulating  to  conducting 
behavior  in  composite  materials.  An  important  feature  of 
the  resulting  model  is  the  indentification  of  quantitative 
structural  parameters  associated  with  the  state  of  orienta¬ 
tion  ( f ^  and/or  f^)  and  filler  shape  as  specified  by  the 
aspect  ratio  {particulate:  a  =  1;  platelet:  a  <  1; 
fiber:  a  >  1).  The  incorporation  of  these  important  struc¬ 

tural  features  necessarily  introduced  computational  complex- 
itites  into  the  model  relationship.  An  interactive  computer 
program  is  described  in  the  appendix  to  facilitate  the  use 
of  the  model  relationships  as  a  screening  tool  for  the 
selection  of  candidate  material  systems. 

For  the  important  special  case  of  relatively  long 
(a  >  20)  highly  conducting  (Kp  >>  Kj^)  fibers  at  volume  frac¬ 
tion  concentrations  greater  than  Vf  =  .1,  the  bulk  conduc¬ 
tivity  can  be  approximated  by  the  following  simplified  form 
of  Equations  (4): 
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Longitudinal:  ”  (l/3)v^  xK°(  1  +  f^)  (1  -  f^) 

Width:  <K2>  =  (l/3)v.  xKpd  -  fp)  (1  "  f^) 

Thickness:  <K2>  =  (l/3)v^  xKp(  1  +  2f^) 

where  is  the  conductivity  of  the  fiber  component.  The 

volume  fraction  concentration,  vf  ,  can  be  obtained  from 

specified  weight  fractions,  Wf  ,  by  Equation  (12)  or  by  the 

calibration  curves  displayed  in  Figure  6.  The  orientation 

parameters,  f  and  f  ,  describe  the  state  of  orientation  of 
a  P 

the  fibers.  Values  of  f^  in  the  vicinity  of  -1/2  indicate 
that  the  fibers  are  in  a  planar  state  of  orientation  with 
the  fibers  confined  to  the  1-2  plane.  The  parameter,  f 

P 

specifies  the  extent  of  orientation  within  the  1-2  plane. 

Values  of  f^  in  the  vicinity  of  1  indicate  fiber  alignment 

along  the  longitudinal  ("1")  axis;  values  of  f  near  0  indi- 

P 

cate  a  planar  random  material.  Three-dimensional  random 

fiber  distributions  are  obtained  v/ith  f  =  f  =0.  The 

a  p 

reduction  factor,  x  >  accounts  for  contact  irregularities 
associated  with  surface  roughness  and/or  surface  coatings 
and  appears  to  be  a  constant  on  the  order  of  1  -  3  x  10~^  . 
The  bulk  resistivity  is  obtained  as  the  reciprocal  of  the 
conductivity;  viz.,  <pi>  =  l/<Ki>  for  i  =  1,2,3. 

The  satisfactory  correlations  obtained  between 
predicted  and  measured  resistivities  lends  confidence  to  the 
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application  of  the  model  relationships  as  a  screening  tool 
for  identifying  candidate  material  systems-  The  rationali¬ 
zation  of  data  scatter  in  terms  of  variations  in  the  state 
of  filler  orientation  emphasizes  the  role  of  microstructure 
in  establishing  transport  properties.  The  dependence  of  the 
bulk  conductivity  on  microstructural  variations  induced  by 
fabrication  procedures  points  out  the  futility  of  attempting 
to  compare  properties  of  unknown  microstructure.  This 
important  conclusion  underscores  the  necessity  of  conducting 
structural  characterizations  in  conjunction  v/ith  traditional 
testing  procedures  for  mechanical  and  transport  properties. 
The  diffraction  techniques  developed  by  McGee  and  McCullough 
[1982]  can  be  used  to  obtain  the  required  structural 


characterizations . 


RE  COMMENDAT IONS 


The  current  model  appears  to  have  captured  the 
major  features  associated  with  electrical  transport  in  com¬ 
posite  materials.  Consequently,  the  model  can  serve  as  a 
preliminary  screening  tool  for  candidate  materials  systems. 
The  relative  succcess  of  these  initial  attempts  to  incor¬ 
porate  structural  parameters  into  relationships  to  predict 
transport  properties  suggests  the  need  for  m.ore  rigorous 
treatments  to  extend  the  model  for  application  to  design 
considerations.  Particular  attention  should  be  given  to 
developing  improved  relationships  to  account  for  the 
influence  of  processing  induced  variations  in  microstructure 
on  the  transport  properties  of  composite  materials. 

The  contact  probability  relationships  used  for 
uniform  spherical  filler  geometries  adequately  predict  the 
critical  concentrations  associated  with  the  insulator- 
conductor  transitions.  These  relationships  should  be 
extended  to  account  for  the  influence  of  particle  size 
distributions  on  the  critical  concentrations.  Literature 
data  suggests  that  the  critical  concentrations  required  to 
achieve  conducting  behavior  can  be  significantly  reduced 
with  the  use  of  broad  distributions  of  particle  size. 
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Although  the  approximations  used  to  estimate  the 
concentration  dependence  of  the  contact  probabilities  for 
fibers  and  platelets  appear  to  capture  the  significant 
trends,  the  treatment  is  an  oversimplification.  Additional 
work  is  required  to  establish  precise  relationships  for 
this  important  structural  feature  of  fiber  and  platelet 
filler  geometries.  Particular  attention  should  be  devoted 
to  incorporating  orientation  effects  into  the  expressions 
for  contact  probabilities. 

The  empirical  reduction  factor,  x  /  used  to 
estimate  the  effective  network  conductance  in  terms  of  the 
filler  conductivity  represents  a  significant  reduction  in 
achievable  conductivities.  A  study  leading  to  the  specifi¬ 
cation  of  the  factors  associated  with  this  reduction  could 
provide  insight  into  surface  and/or  geometric  conditions 
v/hich  could  promote  higher  levels  of  conductance. 

Added  confidence  in  the  model  predictions  could  be 
achieved  through  comparisons  of  predictions  with  an  extended 
data  base.  In  particular,  resistivity  data  for  platelet 
filled  and  multiple  component  systems  (e.g.,  two  or  more 
types  of  filler)  were  not  available  for  test  comparisons 
with  model  predictions.  Consequently,  the  predictions 
obtained  on  such  systems  must  be  viewed  with  some 
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reservation.  It  is  imperative  that  structural  characteriza¬ 
tions  be  conducted  in  concert  with  electrical  and  mechanical 
testing  programs  aimed  at  establishing  an  extended  data 
base.  The  diffraction  technique  developed  by  McGee  and 
McCullough  [1982]  provides  a  method  for  the  characterization 
of  the  state  of  orientation  and  the  situ  aspect  ratio  of 
the  filler.  In  addition,  this  structural  characterization 
technique  can  be  utilized  to  correlate  processing  conditions 
to  the  level  of  conduction  through  the  quantitative  deter¬ 
mination  of  these  processing  sensitive  structural  parameters. 
It  is  reasonable  to  anticipate  that  the  relationship  of  pro¬ 
cessing  conditions  to  material  behavior  can  be  established 
by  linking  these  structure  dependent  models  with  analyses 
designed  to  trace  the  development  of  microstructure  with 
variations  in  processing  conditions. 
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APPENDIX  A 

USER'S  MANUAL  FOR  MICROMECHANICAL  MODELS 
COMPOSITE  MATERIALS 
ELECTRICAL  AND  THERMAL 


RESISTIVITY  AND  CONDUCTIVITY 
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The  following  user's  guide  is  written  to  illustrate 
the  operational  aspects  of  the  computer  program  which  con¬ 
tains  the  theoretical  models  for  the  prediction  of  the  elec¬ 
trical  properties  of  composite  materials.  The  program  is 


I 


written  in  ANSI  standard  FORTRAN-4  with  the  ability  to 
accept  data  as  free  field  input.  This  option  allows  input 
to  the  program  to  occur  anywhere  on  the  input  field  without 
reference  to  column  position.  A  computer  graphics  capabil¬ 
ity  is  implemented  in  the  program  for  the  convenient  display 
of  material  properties.  The  graphics  FORTRAN  call  state¬ 


ments  are  standard  CALCOMP®  statements  and  may  be  imple¬ 
mented  for  interactive  graphics  terminals  or  remote  graphics 
output  devices. 


Two  example  executions  of  the  computer  program  are 
presented.  To  indicate  the  expected  input  and  the 
programmed  flexibilities  inherent  in  the  code,  detailed 
descriptions  written  in  italics  follow  each  interaction  of 
the  first  example.  The  user  input  is  indicated  with  an 
arrow,  . 


For  the  first  example,  we  assume  a  two  phase  system 
and  calculate  the  electrical  properties  for  a  range  of  fiber 
volume  fractions.  The  program  execution  is  given  as 


follows . 
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^run  navy/maincd 
# RUNNING  1781 

The  above  expression  is  the  expected  command  for  a  Burroughs 
B770Q  mainframe.  •  Please  respond  with  an  appropriate  command 
for  each  individual  system. 

CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAL  MODEL  FOR  CALCULATING; 

COMPOSITE  MATERIAL  ELECTRICAL  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIVITIES 

VERSION:  051182 

The  program  responds  with  the  above  heading,  which  indicates 
the  type  of  transport  properties  it  will  calculate  and  its 
version  The  version  number  indicates  that  this  copy  of 

the  program  was  compiled  on  May  11,  1982. 

ENTER>  CALCULATION  OPTION 

(1)  CALCULATE  ELECTRICAL  PROPERTIES 

(2)  CALCULATE  THERMAL  PROPERTIES 

-^1 


The  program  includes  the  abilities  to  calculate  material  proper¬ 
ties  for  electrical  and  thermal  transport  phenomenon.  The  user 
selects  from  the  menu  and  responds.  If  a  number  that  is  not 
contained  in  the  menu  is  entered,  the  above  menu  list  is 
redisplayed.  This  is  the  case  for  all  menu  input  for  this 
program.  The  actual  transport  phenomenon  nvdeled  adjusts  the 
calibration  constants  a  and  8 

ENTER>  INPUT  OPTION 

(1)  INPUT  RESISTIVITY  DATA 

(2)  INPUT  CONDUCTIVITY  DATA 

^1 

For  either  case  of  electrical  or  thermal  transport  properties, 
the  program  will  allow  input  of  either  resistivity  or  conduc¬ 
tivity  constants. 

DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 

The  user  responds  with  either  "NO"  (or  carriage  return)  or  "YES” 
(or  "Y").  Nhen  a  filler  phase  is  requested,  the  program  con¬ 
structs  a  surrogate  matrix  of  resin  and  filler.  At  present,  the 
aspect  ratio  of  the  filler  phase  is  assumed  to  be  unity. 
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TWO  PHASE  SYSTEM  ASSUr4ED:  RESIlY/FILLER 

ENTER>  RESIN  RESISTIVITY 
^4.5E16 


ENTER>  FIBER  RESISTIVITIES 

{ LONGITUDINAL , TRANSVERSE ) 
^1. 5E-3,1.5E-3 


The  ussT  is  prorrrptsd  to  enter  the  isotropic  resin  resistivities 
and  orthotropia  fiber  resistivities.  The  longitudinal  direction 
corresponds  to  the  axial  direction  of  the  fiber,  and  the  trans¬ 
verse  corresponds  to  any  direction  perpe'ndicuZavr  to  the  fiber 


axus. 


ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE,  END  VALUE,  INCREMENT 
->0 . 0 , 1 . 0 , 0 . 1 

Only  for  the  case  of  two  p'nase  composites  will  this  be  displayed. 
This  allows  the  user  to  observe  the  vajn-ance  of  material  proper¬ 
ties  as  the  amount  of  fiber  loading  is  varied..  For  the  oa.se  of  a 
three  phase  corrposite,  this  report  is  not  given.  The  user  may 
select  one  set  of  values  for  resin,  filler  a.nd  fiber  volume 
fraction. 

ENTER>  ASPECT  RATIO  OF  FIBERS 

^20. 


The  aspect  ratio  is  defined  as  the  length  to  diameter  ratio  of 
the  fiber.  For  aspect  ratios  approaching  10^  the  fiber  system 
resembles  continuous  fiber  composites .  Aspect  ratio  equal  to 
one  indicates  spheres.  The  significance  of  c,speat  ratios  less 
than  one  controls  the  program  to  model  platelets.  The  primary 
direction  of  a  fiber  with  aspect  ratio  greater  than  one  will  be 
the  "l"-axis  direction.  For  a  vlatelet  ghase  case,  the  uerven- 
dtculccr  normal  to  the  platelet  surface  will  correspond  with  the 
"1"  axis  direction. 

ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE:  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC:  THREE-DIMENSIONAL 


The  fibers  for  the  in-plane  averaging  will  remain  in  the 
plane  of  the  rrricro-lami'nate  system. 


"1-2" 
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ENTER>  HERMANS  ORIENTATION  PARAMETER 
->.4 

For  the  planar  ty-pe  of  averaging j,  the  variance  of  the  degree  of 
alignment  is  controlled  by  the  Herman's  orientation  parameter 
with  a  range  of  zero  to  one.  The  case  of  "perfectly  random"  is 
represented  by  zero  and  "perfectly  aligned"  is  represented  by 
one.  Most  engineering  materials  will  have  a  value  between  these 
ideal  numbers. 

DEFAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<DATA  CHECK > 

CALIBRATION  PARAMETERS 
ALPHA  =  0.4100 

BETA  =  0.2900 

DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 

These  calibration  parameters  whould  not  need  to  be  changed.  If 
it  is  necessary,  the  program  will  prompt  the  user  by  responding 
with  "YES"  (or  "Y") . 

-^YES  DO  YOU  WISH  A  TABLE  OF  DATA? 


■^NO  DO  YOU  WISH  A  COMPUTER  PLOT? 

These  questions  allow  the  user  to  control  the  flaw  of  output. 
ELECTRICAL 

RESISTIVITIES  OF  COMPOSITE  MATERIAL  :  PLANAR  CASE 


VF 

F 

Kll 

K22 

K33 

0.1000 

0.4000 

3. 2451E+15 

6.7675E+15 

3.6401E+16 

0.2000 

0.4000 

1.0714E-02 

2. 5000E-02 

2.8800E+16 

0.3000 

0.4000 

7.1429E-03 

1.6667E-02 

2.2050E+16 

0.4000 

0.4000 

5.3571E-03 

1. 2500E-02 

1.6200E+16 

0.5000 

0.4000 

4.2857E-03 

l.OOOOE-02 

1.1250E+16 

0.6000 

0.4000 

3.5714E-03 

8. 3333E-02 

7.2000E+15 

0.7000 

0.4000 

3.0612E-03 

7.1429E-03 

4.0500E+15 

0.8000 

0.4000 

2.6786E-03 

6 . 2500E-03 

1. 8000E+15 

0.9000 

0.4000 

2. 3810E-03 

5.5556E-03 

4.5000E+14 

1.0000 

0.4000 

1.5000E-03 

1. 5001E-03 

1. 5001E-03 

The  electrical  resistivities  of  the  composite  material  are 
presented  in  tabular  form.  The  "VF"  column  signifies  the  fiber 
volume  fraction  of  the  fiber  phase.  The  "F"  is  the  Herman's 
orientation  parameter  which  indicates  the  degress  of  randomness 


98 


in  the  composite.  The  "Kjj",  ”-<22"  "-ss"  '^^''■cintities 

vepTesent  the  electrioa.t  resistivities  for  the  longitudi'nal , 
transverse  and  perpendicular  directions^  respectively . 

ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAJyiE  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 

The  user  may  d.irect  the  flew  of  t'^ae  program  for  additional 
calculations  or  to  erud  the  program  execution. 


The  second  example  illustrates  the  calculation  of 
electrical  resistivity  for  an  axially  symmetric  composite. 
The  constitutive  properties  are  assumed  to  be  identical  to 
the  properties  in  the  first  example.  The  program  execution 
is  given  as  follows. 


->run  navy/mained 
#RUNNING  1809 

CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAL  MODEL  FOR  CALCULATING: 
CO.^IPOSITE  iMATERIAL  ELECTRICAL  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIVITIES 

VERSION:  051182 

ENTER>  CALCULATION  OPTION 

(1)  CALCULATE  ELECTRICAL  PROPERTIES 

(2)  CALCULATE  THERMAL  PROPERTIES 

^1 

ENTER>  INPUT  OPTION 

(1)  INPUT  RESISTIVITY  DATA 

(2)  INPUT  CONDUCTIVITY  DATA 

^1 
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-^NO  DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 

TWO  PHASE  SYSTEM  ASSUMED:  RESIN/FILLER 

ENTER>  RESIN  RESISTIVITY 
->4.5E16 

ENTER>  FIBER  RESISTIVITIES 

( LONG I TU  DI NAL , T  RANS  VE  RSE ) 

-^1.5E-3,1.5E-3 

ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE,  END  VALUE,  INCREMENT 

^0.0, 1.0, 0.1 

ENTER>  ASPECT  RATIO  OF  FIBERS 
->20. 

ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE;  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC;  THREE-DIMENSIONAL 

^2 

The  user  has  selected  the  type  of  averaging  for  an  axially 
syrmetria  composite.  The  following  calculations  indicate  this 
type  of  approximation. 

ENTER>  HERMANS  ORIENTATION  PARAMETER 
^0.4 


The  interpretation  of  the  Herman's  orientation  parameter  is 
different  for  the  axially  symmetric  case  than  for  the  planar 
case.  The  Herman’s  orientation  pararr.eter  may  vary  from  -O.S  to 
1.0.  The  value  of  1.0  indicates  a  perfectly  aligned  composite. 

A  three- dimensio'nal  random  composite  is  modeled  with  a  0.0  value 
and  a  transversely  isotropic  composite  with  a  -O.S  value. 

DEFAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<DATA  CHECK> 

CALIBRATION  PARAMETERS 
ALPHA  =  0.4100 

BETA  =  0.2900 

->NO  DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 

->YES  DO  YOU  WISH  A  TABLE  OF  DATA? 


^NO  DO  YOU  WISH  A  COMPUTER  PLOT? 
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ELECTRICAL 

RESISTIVITIES  OF  COMPOSITE  MATERIAL  :  AXIAL  CASE 


VF 

F 

Kll 

K22 

K33 

0.1000 

0.4000 

3. 7305E+15 

9 . 2879E+15 

9.2879E+15 

0 . 2000 

0.4000 

1.2500E-02 

3. 7500E-02 

3 . 7500E-02 

0.3000 

0.4000 

8.3333E-03 

2. 5000E-02 

2.5000E-02 

0.4000 

0.4000 

6.2500E-03 

1.8750E-02 

1.S750E-02 

0.5000 

0.4000 

5. OOOOE-03 

1. 5000E-02 

1.5000E-02 

0.6000 

0.4000 

4.1667E-03 

1. 2500E-02 

1. 2500E-02 

0. 7000 

0.4000 

3.5714E-03 

1. 0714E-02 

1.0714E-02 

0.8000 

0.4000 

3.1250E-03 

9. 3750E-03 

9.3750E-03 

0.9000 

0.4000 

2. 7778E-03 

8. 3333E-03 

8.3333E-03 

1.0000 

0.4000 

1.5001E-03 

1.5001E-03 

1.5001E-03 

The  tabulav  data  pTesented  represents  an  axially  symmetric 
composite.  Note  that  the  K22  electrical  resistivities 

are  equivalent  for  this  material  geometry. 

ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 

The  user  may  direct  the  flow  of  the  program  for  additional 
calculations  or  to  end  the  program  execution. 


APPENDIX  B 

SAMPLE  EXECUTIONS  OF 
MICROMECHANICS  MODELLING 
FOR  ELECTRICAL  RESISTIVITY 
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CONTENTS 


1)  Graphite  Fibers  and  PPS  Resin . 103 

2)  Graphite  Fibers  -  Round  Copper  Particles  and  PPS 

Resin . 115 

3)  Graphite  Fibers  -  Copper  Flakes  and  PPS  Resin  .  .  .  121 

4)  Graphite  Fibers  -  Copper  Fibers  and  PPS  Resin  ...  127 

5)  Nickel  Coated  Graphite  Fibers  and  PPS  Resin  ....  131 


6)  Copper  Fibers  and  PPS  Resin . 141 

7)  Copper  Flakes  and  PPS  Resin . 151 
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1 )  Graphite  Fibers  and  PPS  Resin 

Fiber 

Resistivity  (ohm-cm)  1.5x10“^ 

Volume  Fraction  0.20 

Aspect  Ratio  50 

Orientation  Parameter, 
fp  ,  (planar)  0.5 

Electrical  Resistivities 
Ki  =  1.00x10“^ 

K2  =  3.00x10“^ 

K3  =  2.88x10^® 

The  sample  execution  follows  in  Table  1 .  Another 
execution  is  given  in  Table  2  for  a  scan  over  volume 
fractions  of  the  fiber.  Sample  computer  plots  are  given 
in  Figures  1,  2,  and  3  for  ,  K2  ,  and  K3  , 
respectively. 

Note:  The  reduction  factor,  x  >  been  taken  as  unity. 

The  influence  of  geometrical  irregularities  under 
surface  coatings  may  be  taken  into  account  by 
replacing  the  fiber  resistivity  by  0.45  ohm-cm 
[=  (1.5x10-3)  X  (3x102)]. 


Resin 

4.5x10^^ 

0.80 
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Table  1  (continued) 


DO  YOU  WISH  TO  CHANGE  THESE  DAI 


DO  YOU  WISH  A  TABLE  OF  DATA? 


DO  YOU  WISH  A  COMPUTER  PLOT? 


ELECTRICAL 

RESITIOITIES  OF  COMPOSITE  MATERIAL 
vr  F  Kll 

0»2000  0.5000  l.OOOOE-02 

enter:;-'  repeat  option 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOf 
<3>  START  O'v'ER 

(:4>  END  PROGRAM 


5 
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Table  2. 

Samele  Execution  of  Micrornechanic 
H  o  d  e  1.  f  o  r  G  r  3  h  s  o  f  Electrics  1  R  e  s  i  s  t.  i 
f  o  r  G  r  B  P  h  i  t  e  Fibers  e  n  c!  F'  F'  S  Fv  e  s  i  n 


CENTER  FOR  COMPOSITE  MATERIAI..S 
COMPUTER  ^SuFTUARF  SYSTEMS 

MICROMECHANICAI..  MODEL  FOR  CALCULATING  t 
COMPOSITE  MATERIAL  ELECTICAF  AND  THERMAI. 
RFSISTIUITIES  AND  COMDUCTI VITES 

version:  051182 


ENTER>  CALCULATION  OPTION 

Cl)  CALCULATE  ELECTRICAL  PROPERTIES 
(2)  CALCULATE  THERMAL.  PROPERTIES 


ENTER>  INPUT  OPTION 

(1)  INPUT  RESISTIVITY  DATA 

(2)  INPUT  CONDUCTIVITY  DATA 


DO -YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 


TWO  PH  A 

F:  NTEFv‘> 
<5E"M6 

ENTER> 

<0015?  0 

ENTER> 

<•  0 .?  1 1 0  ? 
ENTER> 


SE  SYSTEM  ASSUMED  1  RESIN/FILLER 
RESIN  RESISTIVITY 


FIBER  RESISTIVITIES 
<  LONGITUDINAL  ?  TRANSVERSE ) 

*0015 

FIBER  VOLUME  FRACTIONS  TO  SCAN 
START  VALUE.?  END  VALUE?  INCREMENT 


ASPECT  RATIO  OF  FIBERS 


ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE  5  TWO-DIMENSIONAL 
( 2  )  AXIAL  S  ¥  M  M  E  T  R I C  ?  T  Fi  R  F  E  -  D I M  E  N  S 1 0  N  A 1.. 
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Table  2  (continued) 


ENTER:>  HFRHA^4S  ORIENTATION  PARAMETER 
><0  '  . 

IiEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<riATA  CHECK> 

CAl-IBRATION  PARAMETERS 
AI..PHA  =  0*4100 

BETA  =  0*2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 


DO  YOU  WISH  A  TABLE  OF  DATA? 

Y 

DO  YOU  UI3H  A  COMPUTER  PLOT? 

Y 


ENTER’--'  GRAF'HIC  OUTF’UT  DAiA 

SELECT  COORDINATE  DIRECTION 
(1)  Kll 
<2)  K22 
<3)  K33 
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Table  2  (continued) 

ELECTF-ICAl... 


ESITIVITIES 

OF  COMPOSITE  HATERIAI..  1 

PLANAR  CASE 

VF 

r~ 

r 

Kll 

K22 

K2>3 

0.0200 

0.5000 

.  3  u  2  /  ii.  T  jL  5 

9.7008E-i-15 

4.3216E+16 

0 .0400 

0  +  5000 

1 .90S5E-H5 

5.2423E-!-15 

4. 1467E-f-16 

\}  +  0  w  0  U 

0,5000 

1 .2389E+15 

3.4985E-1-15 

3  4  7  7  5  5  E  -M  6 

0.0800 

0.5000 

8.'9701E+14 

2.5700E+15 

3.8079E+16 

0 ♦ 1000 

0.5000 

1 . 3  9  2  0 E  4*  0  3 

4. 1759E+03 

3.6450E+16 

0.1200 

0 . 5000 

1 ,6667E-02 

5.0000E-02 

3.4848E-M6 

0.1400 

0.5000 

1 .4236E-02 

4.2857E-02 

3.3282E+16 

0.1600 

0.5000 

1 .2500E-02 

3.7500E~02 

3. 1752E-!-16 

0.1800 

0  4  5000 

l.llllE-02 

3.3333E-02 

3.0258E+16 

0.2000 

0 . 5000 

1 ,0000E-02 

3.0000E-02 

2.8300E+16 

0.2200 

0 . 5000 

9.0909E-03 

2.7273E-02 

2.7378E-M6 

0.2400 

0  ♦  5000 

8.3333E-03 

2 .5000E-02 

2.5992E+16 

0  A  2600 

0 , 5000 

7.6923E-03 

2.3077E-02 

2.4642E+16 

0.2800 

0.5000 

7, 1429E-03 

2.1429E-02 

2.332SEfl6 

0  <-  3000 

0 . 5000 

6  ♦  6667’E-~03 

2  4  OOOOE-02 

2.2050eil6 

0.3200 

0.5000 

6.2500E-03 

1 .8750E-02 

2 . 08088416 

0 .340  0 

0.5000 

5.8824E-03 

1 ,7647E-02 

1 .9602E-1-16 

0 . 3600 

0  ♦  5^000 

5,5556E-03 

1 4  6  6  6  7  E  “  0  2 

1 .34328416 

0 . 3800 

0.5000 

5.2632E-03 

1 ,5789E-02 

1 .7298E416 

0. 4000 

0.5000 

5.0000E-03 

1 .5000E-*02 

1 ,6200E416 

0<  4200 

0,5000 

4.7619E-03 

1 .4236E-02 

1 .5138E416 

0  <  4400 

0.5000 

4,5455E-03 

1 .3636E-02 

1 .41128416 

0 . 4600 

0  4  5000 

4.347SE-03 

1 . 3043E-02 

1 .31228416 

0.4800 

0.5000 

4,1667E-03 

1 .2500E-02 

1 .2168E416 

0 . 5000 

0.5000 

4.0000E-03 

1 . 2000E-02 

1 . 1250E416 

0.5200 

0,5000 

3.S462E-03 

1 . 1538E-02 

1 .0368E416 

0.5400 

0.5000 

3,7037E-03 

1 , llllE-02 

9,5220E415 

0  <  5600 

0.5000 

3 ♦ 5714E-03 

1 .0714E-02 

8.7120E415 

0.5800 

0,5000 

3.4483E-03 

1 .0345E-02 

7.9380E415 

0 . 6<  0  0  0 

0.5000 

3 . 3333E-03 

1 . OOOOE-02 

7.2000E415 

0 » 6200 

0.5000 

3.2258E-03 

9.6774E-03 

6.49S0E415 

0.6400 

0.5000 

3 . 1250E-03 

9.3750E-03 

5'.8320E415 

0 . 6-600 

0 . 5000 

3.0303E-03 

9.0909E— 03 

5.2020E415 

0 . 6800 

0.5000 

2.9412E-03 

O  .  Q  L  “*  U 

4.6080E415 

0.7000 

0.5000 

2.S571E~03 

S.5714E~03 

4.0500E415 

0 . 7200 

0,5000 

2.7778E-03 

3  4  3333E-03 

3.5280E415 

0.7400 

0.5000 

2.7027E-03 

3. 1031E-03 

3.0420E415 

0.7600 

0.5000 

2.6316E-03 

7.S947E~03 

2.5920E415 

0 . 7300 

0.5000 

2.5641E-03 

7.6923E-03 

2.1780E415 

0.3000 

0.5000 

2. 500 OE -03 

7  4  5000E-**03 

1 .8000E415 

0.3200 

0.5000 

2.4390E-03 

7  4  3171£-*03 

1 .4580E415 

0 .34  0  0 

0.5000 

2.3810E-03 

7. 1429E-03 

1-.  1520E415 

0.8600 

0.5000 

2.3256E-03 

6.9767E-03 

3.8200E414 

0.8800 

0 . 5000 

2,2727E-03 

6  4  S 1 3  2^  E 0  3 

6. 4S00E414 

0<9000 

0,5000 

2.2222E-03 

646667E-03 

4.5000E414 

0.9200 

0,5000 

2, 1739E-03 

6.5217E-‘03 

2,8800E414 

0 .9400 

0.5000 

2. 1277E-03 

643330E-03 

1 .6200E4i4 

0c960p 

0.5000 

2. 03.33  E- 03 

642500E-03 

7.2000E413 

0.9800 

0.5000 

2 . 0403E-03 

6. 1224E-03 

1 .S000E413 

1 .0000 

0  4  5000 

1 . 5000E-03 

1 45001E-03 

1 , 5001E-03 
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Table  2  (continued) 


RHADY  FOR  PLOT? 

Y 

ENTER>  OERTICLE  AXIS  MODE 

<1)  PLOT  LINEAR  SCAI.E  FOR  LOniO(Kii) 

(2)  PLOT  LOG  CYCLE  SCALE  FOR  L0610<Kii) 

1 

ENTER>  GRAPHICS  OUTPUT  DEVICE 
(.1.)  TERMINAL 
<2>  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

<  3 )  REPEAT  SAME  PLOT 

<2)  ENTER  NFU  MATERIAL  GEOMETRY 

<3)  START  OVER 

<4)  END  PROGRAM 

n 

ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VAl.UEf  END  VALUE?  INCREMENT 

0.0?  1,0?  0»02 

ENTER>  ASPECT  RATIO  OF  FIBERS 
UO.O 

ENTER>  TYPE  OF  AVERAGING 

(1>  IN-PLANE?  TWO-DIMENSTONAI. 

(2)  AXIAL.  SYMMETRIC?  THRFE-DIMENSIONAL 

1 


ENTER>  HERMANS  ORIENTATION  PARAMETER 
O.U 

DEFLAULT  ELECTRICAL.  CALIBRATION  PARAMETERS 

<DATA  CHECK> 

CALIBRATION  PARAMETERS 
AL.PHA  =0.4100 
BETA  =  0.2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES'? 
N 

DO  YOU  WISH  A  TABLE  OF  DATA? 

N 


Table  2  (continued) 

DO  YOU  WISH  A  COMPUTER  PLOT' 


ENTER>  GRAPHIC  OUTPUT  DATA 

B ELECT  COORDINATE  DIRECTION 
(1.)  Kll 

(2)  K22 

( 3 )  K33 


READY  FOR  PLOT? 


PNTER>  UERTICLE  AXIS  MODE 

(1.)  PLOT  LINEAR  SCALE  FuR  LOGIOCK  J.  i) 
(2)  PLOT  LOG  CYCLE  SCAi_E  FOR  LOGIC  (Kii 


ENTER'?  GRAPHICS  OUTPUT  DEC  ICE 

(1)  TERMINAL 

>'2)  FLATBED  PLOTTER 

;i. 

ENTER?  REPEAT  OPTION 

< 1 )  REPEAT  SAME  PLOT 

(2)  ENTER  NEU  MATERIAL  GEOMETRY 
<3>  START  ODER 

<4)  END  PROGRAM 

ENTER?  FIBEiR:  VOLUME  FRACTIONS  TO  SuAN 

START  VAI.IJE?  END  VAi..UE?  INCREMENT 

0 V  0  ?  1 ♦ 0 ’  0  »  02 

ENTER?  ASPECT  RATIO  OF  FIBERS 
00*0 

ENTER?  TYPE  OF  AVERAGING 

(1  )  IN~PLANE?  TWO-DIMENSIONAl. 

(2)  AXIAL  SYMME'IRICf  THREE-DIMENSIONA! 
I 

ENTER?  HERMANS  ORIENTATION  PARAMETER 

0  4  5 

DEFLAULT  ELECTRICAI-  CALIBRATION  PARAMETERS 

ADATA  CHECK? 

CAI..IBRATION  PARAMETERS 
lA  1 ..  F'  H  A  0*41  'J  u 

DF.  )  A  “  0  ♦  2Y0u 

DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 

M 

DO  YOU  WISH  A  TABLE  OF  DHTHr 


DO  YOU  WISH  A  COMF'UTER  PLOT? 


Table  2  (continued) 


ENTER>  GRAPHIC  OUTPUT  DATA 

SF.I.ECT  COORDINATE  DIRECTION 
<1 >  Kll 
(?)  K22 
<  3 )  K33 


READY  FOR  PLOT? 


li-NTER^-  VERTICLE  AXIS  MODE 

PLOT  LINEAR  SCAI  E  FOR  LOGlO(Kii) 

<2>  PLOT  LOG  CYCLE  SCALE  FOR  LOGlOCKii) 


ENTER>  GRAPHICS  OUTPUT  DEVICE 
(1)  TERMINAL 
<2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

<1)  REPEAT  SAME  PLOT 

<2)  ENTER  NEW  MATERIAL.  GEOMETRY 

<3)  START  OVER 

<4)  END  PROGRAM 


Figure  1.  Electrical  resistivity  for  graphite/PPS  composite,  longitudinal  direction. 


Figure  2.  Electrical  resistivity  for  graphite/PPS  composite,  transverse  direction. 
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ber  Volume  Fraction^  Vf 


Figure  3.  Electrical  resistivity  for  graphite/PPS  composite,  perpendicular  direction. 
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2 )  Graphite  Fibers  -  Round  Copper  Particles  and  PPS  Resin 


Fl  " 

Fiber 

F2  ^ 

Copper 

Particles 

M 

Resin 

Resistivity  (ohm- 

cm) 

1.5x10"^ 

l.TOxlO”® 

4.5x10^® 

Volume  Fraction 

0.30 

0.20 

0.50 

Aspect  Ratio 

50 

1.0 

— 

Orientation  Parameter, 

fa  ,  (axial) 

0.5 

0.0 

— 

Calculations 

Step  I:  Apparent  volume 
fibers  and  PPS 

fractions 

resin. 

for  System  I 

graphite 

V*  - 

0.30 

=  0.375 

^1  1 

1.0  -  0.2 

V*  -  _ 

0.50 

=  0.625 

''f  1 

mi 

-''^2 

1.0  -  0.2 

Resistivities  for 

System 

I  (  see  Table  3  ) 

-3 

6.00x10 

ohm-cm 

-2 

2.40x10 

ohm-cm 

-2 

2.40x10 

ohm-cm 

Note:  The  influence  of  geometrical  irregularities  and/or 

surface  coatings  may  be  taken  into  account  by 
replacing  the  filler  resistivity  by  the  effective 
resistivity. 

Effective  resistivity  of  graphite  fiber  =  0.45  ohm-cm. 

^Effective  resistivity  of  copper  particles  =  1.7xl0~3  ohm-cm. 
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Step  II:  Apparent  volume  fractions  for  System  II 
Round  Copper  Particles  and  PPS  Resin 


"2 

1  -  V, 


0.20 

1.0  -  6.30 


=  0.286 


1  -  V, 


0.50 

1.0  -  0.30 


=  0.714 


Resistivities  for  System  II  (see  Table  4) 


1  16 
Pjj  =  1.784x10  ohm-cm 


Pjj  =  1.784x10^'^  ohra-cm 


3  14 

Pj^  =  1.784x10  ohm-cm 


Step  III:  Calculate  properties  of  three-phase  system 

pi  =  (1  -  Vi  )pi  t  (1  -  V'  -  Vj  P„ 

2  1  m 


l/P  =  (1  -  0.286)  (- 


)  +  (1  -  0.375) 


6  .00x10 


-)  -  (0.50)  (- 


1 .784x10 


4.5x10' 


p^  =  8.403x10”^  ohm-cm 


1/P^  =  (1  -  0.286)  ( - - - 

2 .40x10' 


)  +  (1  -  0.375) 


-)  -  (0.50)  (- 


1.784x10 


4.5x10' 


2  3  -2 

p  =  P  =  3.36x10  ohm-cm 
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Table  3. 

Calculations  for  SysteiTi  I  for 
Graphite  Fibers  and  PPS  Resin 


CENTER  FOR  COMPOSITE  nATERIAI..S 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAL  MODEL  FOR  CALCULATING T 
COMPOSITE  MATERIAL  ELECTICAI  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIVITES 

VERSION J  051182 


ENTER>  CAi..CULATION  OPTION 

(1)  CAICULATE  ELECTRICAL  PROPERTIES 

(2)  CAI..CULATE  THERMAL  PROPERTIES 


ENTER>  INPUT  OPTION 

Cl)  INPUT  RESISTIVITY  DATA 
■  (2)  INPUT  CONDUCTIVITY  DATA 


DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 


TWO  PHASE  SYSTEM  ASSUMED!  RESIN/FILLER 

1VNTER>  RESIN  RESISTIVITY 
-^•<5E+16 

enter;;:-  fiber  resistivities 

(LONGITUDINAL  j  TRANSVERSE ) 

0<00i5?  0.0015 

enter:;:-  fiber  volume  fractions  TO  SCAN 

START  VALUE?  END  VAI.UE?  TNCREiJIENT 
0<375?  0.375?  1.0 

enter:;:-  aspect  ratio  of  fibers 

50.0 

enter:;:-  type  of  averaging 

(1)  IN-PLANE?  two-dimensional. 

(2)  AXIAL  SYMMETRIC?  THREE-DIMENSIONAL 


enter:;:-  HERMANS  ORIENTATION  PARAMETER 
0.5 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

-::;data  chfck;;:- 

CALIBRATION  PARAMETERS 
AI..PHA  =  0.4100 

BETA  =  0.2900 
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Table  3  (continued) 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 
N 


DO  YOU  WISH  A  TABLE  OF  DATA?' 


DO  YOU  WISH  A  COMPUTER  PLOT? 


ELECTRICAL 

RFSITIVITIES  OF  COMPOSITE  MATERIAL 


A'XIAL  CASE 


VF 

0<  3750 


F 

0 . 5000 


Kll 

A ♦ OOOOE-03 


K22 

2.4000E-02 


2»4000F-02 


EH"lER>  REPEAT  OPTION 

t;  ;l  )  REPEAT  SAME  PLOT 

(2)  ENTER  NFU  MATER  I  Al.  GEOMETRY 

(3)  START  OVER 
<4>  END  PROGRAM 


119 


Tsble  4v 

Cslculstions  for  System  II  for 
Round  Copper  Particles  end  PPS  Resin 


CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAL  MODEL  FOR  CALCULATINGS 
COMPOSITE  MATERIAL  ELECTICAl..  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIVITES 


VERSIONS  051182 


ENTER> 


1 


CAL 
( 1. ) 
(2) 


CULATION  OPTION 

CALCULATE  ELECTRICAL.  PROPERTIES 
CALCULATE  THERMAL  PROPERTIES 


ENTER'.;-  INPUT  OPTION 

<1)  INPUT  RESISTIVITY  DATA 
<2)  INPUT  CONDUCTIVITY  DATA 

1 

DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 
N 

TWO  PHASE  SYSTEM  ASSUMED  S  RESIN./FILLFR 

1VNTER;:=-  RESIN  RESISTIVITY  ’ 

4,5E+16 

ENTER>  FIBER  RESISTIVITIES 

<  LONG  I  TUDINAL.»  TRANSVERSE) 
lt7E—6>  1*7E— 6 


ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VAt.UE  ?  END  VAL  UE  »  .  INCREMENT 
0.?.Q6f  0*2S6?  1.0 

ENTER>  ASPECT  RATIO  OF  FIBERS 
:i.  .0 

enter;:-  type  of  averaging 

<1)  IN-PLANE5  TWO-DIMENSIONAL 

(?)  AXIAL  SYMMETRIC?  THREE-DIMENSIONAL. 
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Table  4  (continued) 


EN!TER>  HERHANS  ORIENTATION  PARAHFTER 
•0 , 0 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

CriATA  CHFCK> 

CALIBRATION  PARAMETERS 
AI.PHA  =  0.4:1.00 

BETA  =  0.:2900 


DO  YOU  NISH  TO  CHANGE  THESE  VAiJJES? 


DO  YOU  WISH  A  TABLE  OF  DATA? 


DO  YOU  WISH  A  COMPUTER  PLOT? 


ELECTRICAL 

RFSITIOITIES  OF  COMPOSITE  MATERIAL  t  AXIAL  CASE 


OF  F  Kll  K22  K33 


0.2360  0.0000  1.7839E+16  1.7839E+16 

ENTER;>  REPEAT  OPTION 

<1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  ODER 
<4)  END  PROGRAM 
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Graphite  Fibers  -  Copper 

Flakes  and 

PPS  Resin 

Fl  ^ 

F2  ^ 

Copper 

M 

Fiber 

Particles 

Resin 

Resistivity  (ohm-cm) 

1.5x10"3 

1 .70x10”® 

4.5x10^® 

Volume  Fraction 

0.2 

0.5 

0.3 

Aspect  Ratio 

50 

1 .0x10"® 

Orientation  Parameter, 

fa  ,  (axial) 

0.5 

0.0 

Calculations 

Step  I ;  Apparent  volume 
fibers  and  PPS 

fractions 

resin. 

for  System  I 

graphite 

V 

V  = 

''f.  1  -  V- 

0.2 

1.0  -  0.5 

=  0.4 

^  ^2 

V  _ 

fn, 

V-  -  - 

0.3 

=  0.6 

1  -  V. 


1.0  -  0.5 


Resistivities  for  System  I  (see  Table  5) 


1  -3 

Pj  =  5.616x10  ohm-cm 


2  -2 
Pj  =  2.223x10  ohm-cm 


3  -2 

Pj  =  2.223x10  ohm-cm 


Note :  The  influence  of  geometrical  irregularities  and/or 

surface  coatings  may  be  taken  into  account  by 
replacing  the  filler  resistivity  by  the  effective 
resistivity. 


Effective  resistivity  of  graphite  fiber  =  0.45  ohm-cm. 
^Effective  resistivity  of  copper  particles  =  1.7xl0“3  ohm-cm. 
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Step  II:  Apparent  volume  fractions  for  System  II 
Copper  Flakes  and  PPS  Resin 


V 


0.5 


1  -  V, 


1.0  -  0.2 


0.625 


V, 


^f 


m 


0.3 


m. 


1  -  V, 


1.0  -  0.2 


=  0.375 


Resistivities  for  System  II  (see  Table  6) 

1  -6 
Pjj  =  4.080x10  ohm-cm 


2  -4 

Pjj  =  4.080x10  ohm-cm 


Pjj  =  4.080x10  ohm-cm 


Step  III:  Calculate  properties  of  three-phase  system 

1/P^  =  (1  -  0.615)  ( - -  )  +  (1  -  0.4) 

5.616x10" 


(- 


4.080x10 


-6 


)  -  (0.3)  (- 


4.5x10 


16' 


1  —  6 

p  =  6.80x10  ohm-cm 


l/P^  =  (1  -  0.625)  { - - -  )  +  (1  -  0.4) 

2.223x10 


(- 


4.080x10 


1  -)  -  (0.3)  ( - ^,) 


-6‘ 


4.5x10 


16' 


2  3  6 
p  =  P  =  6.80x10  ohm-cm 


123 


Tsble  5» 

Ci.5lculste  Resistivities  for  System  I 
Greehite  Fibers  end  PF'S  Resin ♦ 


CFNTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAI..  MODEL  FOR  CALCULATING: 
COMPOSITE  MATERIAL  ELECTICAI..  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTI VITES 

version:  051182 


ENTER>  CALCULATION  OPTION 

.  (1)  CAI  CLILATE  ELECTRICAL  PROPERTIES 

(2)  CALCULATE  THERMAL  PROPERTIES 


enter;:--  input  OPTION 

<1)  INPUT  RESISTIVITY  DATA 
(2>  INPUT  CONDUCTIVITY  DATA 


DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 
ri 

TWO  PHASE  SYSTEM  ASSUMED:  RFSIN/FILLER  / 

enter::=-  resin  resistivity 

4  <  5 


enter:>  fiber  resistivities 

( LONGITUDINAL » TRANSVERSE ) 

0<0015?  0»0015 

ENTER-:-  fiber  VOLUME  FRACTIONS 'TO  SCAN 

rn'ART  VALUED  END  VALUE?  INCREMENT 
0  <  -4  ?  0  ♦  4  f  1  *  0 

enter;::-  aspect  ratio  of  fibers 

50.0 

enter;::-  type  of  averaging 

<1>  IN-PLANE?  TWO-DIMENSIONAI.. 

(2)  AXIAL  SYMMETRIC?  THREE-DIMENSIONAL 


enter;::-  HERMANS  ORIENTATION  PARAMETER 
0.5 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 
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Table  5  (continued) 


CriATA  CHECK> 

C  A 1 ..  I B  R  A  T  T  0  N  R  A  R  A  M  E  T  E  R  S 
ALPHA  ==  0.4100 

••  BETA  ^  0.2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 


DO  YOU  WISH  A  TABLE  OF  DATA? 


DO  YOU  WISH  A  COMPUTER  PLOT? 


ELECTRICAL 

ESITIVITIES  OF  COMPOSITE  MATERIAL  J  AXIAL  CASE 


VF  F  Kll  K22  K33 

0.4000  0.5000  5.6157E-03  2.2232E-02  2.2232E-02 

ENTER>  REPEAT  OPTION 

( 1 )  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIA!..  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 

/. 
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Table  6 4 

Calculate  -Reisistivities  for  System  I 
Coeeer  Flakes  and  PPS  Resin* 


f 

CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROMECHANICAL  MODEL  FOR  CALCULATING? 
COMPOSITE  MATERIAL  ELECT I CAL  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIUITES 

VERSION?  051182 


ENTER>  CALCULATION  OPTION 

< 1 )  CALCULATE  ELECTRICAL  PROPERTIES 
(2)  CALCULATE  THERMAL  PROPERTIES 

1 

ENTER>  INPUT  OPTION 

<t)  INPUT  RESISTIVITY  DATA 
(2)  INPUT  CONDUCTIVITY  DATA 

N  DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 

TWO  PHASE  SYSTEM  ASSUMED?  RESIN/FILLER 

ENTER>  RESIN  RESISTIVITY 
4,5EM6 

ENTER>  FIBER  RESISTIVITIES 

( LONGITUDINAL p TRANSVERSE) 

1*7E-6p  1*7E-6  -  . 

ENTER>  FIBER  VOLUME  FRACTIONS'  TO  SCAN 

START  VALUE?  END  VALUE?  INCREMENT 
0.625*0*625?  1*0 

ENTER>  ASPECT  RATIO  OF  FIBERS 
l*0E-5  ■ 

ENTER>  TYPE  OF  AVERAGING 

(1)  IH-PLANE?  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC?-  THREE-DIMENSIONAI. 


enter;:--  HERMANS  ORIENTATION  PARAMETER 
0*0- 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 
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Table  6  (continued) 


•C'.DATA  CHECIO 
I"  A  L..  I B  R  A  T 1 0  N  P  A  R  A  M  E  T !”  S 
ALPHA  ™  0.>4:!.00 


N  DO  YOU  WISH  TO  CHANGE  THESE  UALUES? 
Y  DO  YOU  WISH  A  TABLE  OF  DATA? 

N  DO  YOU  WISH  A  COMPUTER  PLOT? 

electrical 

RE3ITIVITIES  OF  COMPOSITE  MATE  I 


F  Kll.  K22  K33 

0.6250  0.0000  4.0000E-06'  4.0800E~06  4.0800E-06 

ENTER!  REPEAT  OPTION 

( 1 )  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 
03)  START  OVER 

^'4)  END  PROGRAM 


aAL 


AXIAL  CASE 
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Fi  ^ 

F2  ^ 

M 

Fiber 

Copper 

Fibers 

Resin 

Resistivity  (ohm-cm) 

1.5x10"3 

1.7x10“® 

4.5x10^® 

Volume  Fraction 

0.2 

0.5 

0.3 

Aspect  Ratio 

50.0 

50.0 

— 

Orientation  Parameter, 
fa  ,  (axial) 

0.5 

0.5 

Calculations 

Step  I:  Apparent  volume  fractions  for  System  I  graphite 
fibers  and  PPS  resin. 


0.4 


0.3 

1.0  -  0.5 


0  .6 


Note :  rTie  influence  of  geometrical  irregularities  and/or 

surface  coatings  may  be  taken  into  account  by 
replacing  the  filler  resistivity  by  the  effective 
resistivity. 

& 

Effective  resistivity  of  graphite  fiber  =  0.45  ohm-cm. 
bEffective  resistivity  of  copper  particles  s  1.7x10-3  ohm-cm. 


step  II:  Apparent  volume  fractions  for  System  II 
Copper  Fibers  and  PPS  Resin 


0.5 

1.0  -  0.2 


0.625 


0.3 

1.0  -  0.2 


0.375 


Resistivities  for  System  II  (see  Table  7) 

1  —6 

Pjj  =  4.080x10  ohm-cm 

2  -5 

Pjj  =  1.632x10  ohm-cm 

3  -5 

Pjj  =  1.632x10  ohm-cm 


Step  III:  Calculate  properties  of  three-phase  system 

1/P^  =  (1  -  0.625)  {- - - - -  )  +  (1  -  0.4) 

5.616x10 


( - - - -  (0.3)  ( - ^^  ) 

4.030x10  4.5x10 


p^  =  7.00x10  ^  ohm-cm 


1/P^  =  (1  -  0.625)  ( - i -  )  +  (1  -  0.4) 

2.223x10" 


( - - - ^)  -  (0.3)  ( - ) 

1.632x10  4.5x10 


p^  =  P^  =  2.7x10  ^  ohm-cm 
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Table  . 

Calculste  Resistivities  for  System  II 
Coeeer  Fibers  and  PF‘S  Resin* 


Cr-.1-4TER  FOR  COMPOSITE  MATERIALS 
CiOMPUTER  SOFTWARE  SYSTEMS 

liJCROMECHANICAL  MODEL  FOR  CAI  CULATINGJ 
COMPOSITE  MATERIAL.  ELECTICAL  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTI VITES 

Vr-RSIOHt  051182 


r£NTER:>  CALCULATION  OPTION 

(1)  CALCULATE  ELECTRICAL  PROPERTIES 

(2)  CAIXULATE  THERMAI..  PROPERTIES 

1 

ENTER>  INPUT  OPTION 

<1)  INPUT  RESISTIVITY  DATA 
(2)  INPUT  CONDUCTIVITY  DATA 

1 


DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 
N 


TWO  PHASE  SYSTEM  ASSUMED**  RESIN/FILLER 

ENTER>  RESIN  RESISTIVITY 
4*5E-}-16 

ENTER>  FIBER  RESISTIVITIES 

( LONGITUDINAL » TRANSVERSE) 

1*7E-S.  1.7E-6 

ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE »  END  VAlXEf  INCREMENT 
0*625  3.  0*625f  1*0 

ENTER>  ASPECT  RATIO  OF  FIBERS 
50*0 


ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE?  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC?  THREE-DIMENSIONAL 


enter:;--  HERMANS  ORIENTATION  PARAMETER 
0  <  5 

DEFLAULT  ELECTRICAL  CAI..IBRATION  PARAMETERS 


Table  7  (continued) 
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CDATA  CHF:CK> 

C  A I  I B  Fi;  A  T 1 0  N  P  A  R  A  M  E  T  E  R  S 
A!.  Pi-1  A  =  0,4100 

SETA  =  0,2900 


SO  YOU  WISH  TO  CHANGE  THESE  OALUES? 

DO  YOU  WISH  A  TABLE  OF  DATA? 

V 

DO  YOU  WISH  A  COMPUTER  PLOT? 
electrical 

case 


K11 


0,6250 


u,5u0o  4,0800E-06  1.6320E-05 


enter;:--  REPEAT  OPTION 

<■  1  )  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL,  GEOMETRY 

(3)  START  OVER 

^  4 )  E  N  l!  F’  R  0  G  F-!  A  H 
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5 )  Nickel  Coated  Graphite  Fibers  and  PPS  Resin 

.  a 

Nickel  Coated  Fiber  Resin 

—  16 

Resistivity  {ohm-cm)  2.5x10  4.5x10 

Volume  Fraction  Variable 

Aspect  Ratio  50.0 

Orientation  Parameter, 

fa  /  (axial)  0.5 

•  Executions  are  given  in  Table  8  for  the  computer 

graphs  depicted  in  Figures  4,  5,  and  6  for  ,  K2  , 
and  Kg  ,  respectively. 


^Effective  resistivity  for  damaged  coatings  (see  Table  VIII 
of  text) . 
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Table  8* 

Calculations  for  Nickel  Coated 
Fibers  and  PF'S  Resin  Composite* 


CENTER  FOR  COHPOSITE  MATERIALS 
COMPUTER  SOFTUARF  SYSTEMS 

M I C  R  0  M  E  C  H  A  N I C  A  L  M  0  D  E  L  FOR  C  A  L  C  U 1.  A  T I N  G  J 
COMPOSITE  MATERIAi.^  ELECTTCAI  AND  THERMAL 
RESISTIVITIES  AND  CONDUCTIVITES 

version:  051132 


ENTER>  CALCULATION  OPTION 

(1)  CAi.CULATE  EI..FCTRICAL  PROPER! I 

(2)  CALCULATE  THERMAL  PROPERTIES 


ENTER!  INPUT  OPTION 

(1>  INPUT  RESISTIVITY  DATA 
(2)  INPUT  CONDUCTIVITY  DATA 


f 


DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE 
N 


TWO  PHASE  SYSTEM  ASSUMED:  RESIN/FILLER 

ENTER!-  RESIN  RESISTIVITY 

*  5E  + 1 6  f 

ENTER!-  FIBER  RESISTIVITIES 

( LONGITUDINAL. » TRANSVERSE  ) 

*0025?  0*0025? 

ENTEFL>  F  I. BER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE?  END  VALUE?  INCREMENT 

*  0  ?  1  *  0  ?  0*02? 


ENTER!-  ASPECT  RATIO  OF  FIBERS 
0*0? 
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Tc±)le  8  (continued) 


ENTEFO  TYPE  OF  AVERAGING 

(1)  IN-PLANE5  TUO-DIiiENSIONAL 

(2)  AXIAL  SYMMETRIC?  T1-IREE--DIMENSI0NAL 


ENTER>  HERMANS  ORIENTATION  PARAMETER 


DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<IiATA  CHECIO 
CALIBRATION  PARAMETERS 
ALPHA  ==  0.4100 

BETA  =  0.2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 


•f  .  DO  YOU  WISH  A  TABLE  OF  DATA? 

Y 

•f  DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 

ENTER>  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 
<1.)  Kll 
<2)  K22 

(3)  K33 
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Table  8  (continued) 
ELECTRICAL 


SITIVIT 

lES  OF  COHPOS 

ITE  MATEFLIAL  J 

PLANAR  CASir, 

yP 

F 

Kll 

K22 

K33 

0  *  0200 

0.5000 

3.8027E115 

9.7008E+15 

4.3216E-fl6 

0  *  0400 

0.5000 

1 .9085Efl3 

5*242SE-M5 

4. i467Efl6 

0,0600 

0.5000 

1 .2339E-H5 

3.4983Et15 

3.9753EH6 

0 . 0800 

0 . 5000 

8. 970 IE +14 

2.5700Efl5 

3.8079Efl6 

0.1000 

0.5000 

1 . 6605E+08 

4.98i4E-l-0S 

3,6450E-f  16 

0*1200 

0 , 5000 

2.7778E-02 

8.3333E-02 

3 .4848Ef 16 

0.1400 

0.5000 

2.3810E-02 

7 . 1429E-02 

3.3282E+16 

0.1600 

0  ♦  5000 

2.0B33E-02 

6.,2500E-02 

3. 1752E+16 

0. 1800 

0  *  5000 

1 .85198-02 

5.5556E-02 

3.0253E'H6 

0.2000 

0.5000 

1 .6667E-02 

5. 00 00 E -02 

2.8800Efl6 

0,2200 

0.5000 

1 ,5152E-02 

4  *  5455E-02 

2,7378£+16 

0.2400 

0  *  5000 

1 .3889E-02 

4.1667E-02 

2.5992E+16 

0.2600 

0  ♦  5000 

1 ,2821E-02 

3.8462E-02 

2.4642E+16 

0 . 2800 

0.5000 

1.1905E-02 

3.5714E-02 

2.332SE-fl6 

0.3000 

0*5000 

l.llllE-02 

3.3333E-02 

2 .2050EH6 

0  *  3200 

0 . 5000 

1 .0417E-02 

3,1250E-02 

2.0308Efl6 

0.3400 

0  *  5000 

9, 8039 E -03 

2.,9412E-02 

1 .9602EM6 

0  *  3  6 0 0 

0.5000 

9.2593E-03 

2.7778E-02 

1 ,8432E+16 

0  *  3800 

0 . 5000 

8.7719E-03 

2.6316E-02 

1 .7298E+ 16 

0  *  4000 

0 ♦ 5000 

3  +  333 3  E  ““  0  >:> 

2.5000E-02 

1 *6200Efl6 

0 ♦ 4200 

0 . 5000 

7  *  9365E-03 

2.3S10E-02 

1 .5138E+16 

0 . 4400 

0.5000 

/  *  O  fcl  sJ 

2.2727E-02 

1 .4112E-fl6 

0*4600 

0 ♦ 5000 

7,2464E~03 

2.1739E-02 

1 .3122E-fl6 

o 

CO 

o 

o 

0*5000 

6 . 9444E-03 

2.0833E-02 

1 .2168E+ 16 

0*5  0  0  0 

0*5000 

6 . 6667E— 03 

2.0000E -02 

1 . 1250E+16 

0.5200 

0*5000 

6.4103E-03 

1 .9231E-02 

1 .0368E+16 

0.5400 

0»5i000 

6, 1728E-03 

1 .8519E-02 

9.5220E+15 

0.5600 

0*5000 

5.9524E-03 

1 .7857E-02 

8.7120E+i5 

0*5800 

0  *  5000 

5.7471E-03 

1 .724  IE -02 

7.9330E+15 

0  *  6000 

0*5000 

5.5556E-03 

.1. .  w  ^5  /  E  —  2 

7.2000Eil5 

0*6200 

0*5000 

5 . 3763E-03 

1 .6129E-02 

6,4980E+15 

0.6400 

0*5000 

5.20B3E-03 

1 .5625E  -02 

5.3320E+i5 

0*6600 

0*5000 

5.0505E-03 

1 .5152E-02 

5,2020Et15 

0 . 6800 

0*5000 

4. 9020 E -03 

1.4706E-02  • 

4,6080E+15 

0.7000 

0*5000 

4,7619E-03 

1 .4286E-02 

4.0500E+15 

0 . 7200 

0*5000 

4.6296E-03 

1 .3389E-02 

3 .5280E+15 

0.7400 

0  *  5000 

4 .5045E-03 

1.3514E.-02 

3.0420E+15 

0  *  7600 

0*5000 

4.3860E-03 

1 .3153E-d2 

2.5920Et15 

0.7800 

0*5000 

4.2735E-03 

1 .2 82 IE -02 

2 , 1780ET 15 

0  *  8000 

0*5000 

4 . 1 667 E- 0-3 

1 .2500E-02 

1 .8000E+15 

0.8200 

0  *  5000 

4 , 0650E-03 

1 ,2195E-02 

1 .4580E+15 

0*3400 

0  *  5000 

3.9633E-03 

.  1.1905E-02 

1 . 1520E+15 

0  *  8600 

0  *  5000 

3.8760E-03 

1 .1628E-02 

3.3200E+14 

0  +  8  £^00 

0*5000 

3.7879E-03 

1 . 1364E-02 

6 . 4  3  0  0  E  + 1  4 

0  *  9  0 0 0 

0*5000 

3 . 7037E-03 

1 . llllE-02 

4 .SOOOEf 14 

0*9200 

0*5000 

3.6232E-03 

1 . 0870E-02 

2,3300Efl4 

0*9400 

0*5000 

3.5461E-03 

1 . 0638E-02 

1 ,6200E+14 

0*9600 

0  *  5000 

3 ,4722E-03 

1 .0417E-02 

7.2000E+13 

0*9300 

0  *  5000 

3.4014E-03 

1 .0204E-02 

1 .8000E+13 

1  A  u  0  0  0 

0  *  5000 

2. 500 OE -03 

2 .5001E-03 

2 .500  IE -03 

Table  8  (continued) 
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+  REZADY  FOR  PLOT? 

Y 

ENTEF'O  UERTICLE  AXIS  MODE 

(1)  PLOT  LINEAR  SCALE  FOR  LOGlO(Kii)  ' 
<2)  PLOT  LOG  CYCLE  SCALE  FOR  LOGIC  (K:ii) 


enter::--  graphics  output  device 

(1)  terminai. 

(2)  FLATBED  PLOTTER 

2» 

enter;::-  repeat  OPTION 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  .GEOMETRY 
(3>  START  OVER 

(4)  END  PROGRAM. 


enter;::-,  fiber  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE »  END  VALUE »  INCREMENT 

0»0?  l»0if  0  *  02  7 

enter;-.:-  aspect  ratio  of.  fibers 

50,07  , 


enter;::-  type  of  averaging 

(1)  IN-PLANE}  TWO-DIMENSIONAL 

(2)  AXIAI-  SYMMETRIC?  THREE-DIMENSIONAL 


enter;::-  HERMANS  ORIENTATION  PARAMETER 
0.57 


DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

-::;data  check;::- 

CALIBRATION  PARAMETERS 
ALPHA  =  0.4100 

BETA  =  0.2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 


+  DO  YOU  WISH  A  TABLE  OF  DATA? 
N 

•f  DO  YOU  WISH  A  COMPUTER  PLOT? 
Y 


enter;::-  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 

(1)  Kll 

(2)  K22 
<  3 )  K33 


Table  8  (continued) 


READY  F(3R  PLOT'? 
Y 

ENTER 


1 

ENTER>  GRAPHICS  OUTPUT  DEVICE 

(1)  TERMINAI 
<2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

<1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEU  MATERIAL  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 


ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE?  END  VALUE?  INCREMENT 
0»0?  :i.  »0?  0  +  02? 

ENTER>  ASPECT  RATIO  OF  FIBERS 
5  0  •>  0  ? 

ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE?  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC?  THREE-D IMENS I ONAL 

2  ? 

ENTER>  HERMANS  ORIENTATION  PARAMETER 

0*5? 

DEFLALILT  ELECTRICAL  CALIBRATION  PARAMETERS 

CDATA  CHECIO 
CALIBRATION  PARAMETERS 
ALPHA  ==  0,4100 

BETA  =  0*2900 


•f  DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 
N 

■}■  DO  YOU  WISH  A  TABLE  OF  DATA? 

N 

T  DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 

ENTER>  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 

(1)  Kll 

(2)  K22 

(3)  K33 
.5  ? 


VERTICLE  AXIS  MODE 

<1)  PLOT  LINEAR  SCALE  FOR  LOGlO(Kii) 
(2)  PLOT  LOG  CYCLE  SCAI.E  FOR  LOG:LO(Ki:L 
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Table  8  (continued) 


READY  FOR  PLOT? 
Y 


EHTER>  VERTICLE  AXIS  MODE 

(1)  PLOT  LINEAR  SCALE  FOR  LOGIO(Kii) 

(2)  PLOT  LOG  CYCLE  SCALE  FOR  L0G10<Kii 

1  r 

ENTER>  GRAPHICS  OUTPUT  DEVICE 

(1)  TERMINAi. 

(2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAME  PLOT 

<2>  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 
4f 


Figure  4.  Electrical  resistivity  for  nickel  coated  graphite  fibers  and  PPS  resin, 
longitudinal  direction. 
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Figure  5.  Electrical  resistivity  for  nickel  coated  graphite  fibers  and  PPS  resin 
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Figure  6.  Electrical  resistivity  for  nickel  coated  graphite  fibers 


.00  0.20  0.40  0.60  0.80  1.00 

Fi  b  er  Volume  Fraction^  Vf 
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6 )  Copper  Fibers  and  PPS  Resin 


Resistivity  (ohm-cm) 

Volume  Fraction 

Aspect  Ratio 

Orientation  Parameter, 
fa  ,  (axial) 


Copper  Fibers^ 
1.70x10”® 
Variable 
50.0 


0.5 


Resin 

4.5x10^® 


Executions  are  given  in  Table  9  for  the  computer 
graphs  depicted  in  Figures  7,  8,  and  9  for  ,  K2  , 
and  K3  ,  respectively. 


Note:  The  influence  of  geometrical  irregularities  and/or 

surface  coatings  may  be  taken  into  account  for 
replacing  the  filler  resistivity  by  the  effective 
resistivity. 

^Effective  resistivity  of  copper  2  1.70xl0”3  ohm-cm. 
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Tsbl  (•?  9* 

C  s  1  c  u  1. 3 1  i  o  Tf  s  r  o  r  R  e  s  i  s  t  i  v  i  t  :i  e  •=• 
■for  Copper  Fibers  srid  F’f''S  Resi rit 


CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

MICROHECHANTCAL  MODEL  FOR  CALCULATING: 
COMPOSITE  MATERIAI..  ELECTICAL  AND  THERMAL 
RFSiSTI'v'ITIES  AND  CONDUCTIVITES 


version:  00:11. 


ENTER!  CALCULATION  OPTION 

(1>  CAi..Cl.lLATE  ELECTRICAL  PROPERTIES 
<2)  CA1..CULATE  THERMAL  PROPERTIES 


ENTER!-  INPUT  OPTION 

(1)  INPUT  RESISTIVITY  DATA 

(2)  INPUT  CONDUCTIVITY  DATA 


DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE” 


TWO  PHASE  SYSTEM  ASSUMED::  RESIN./I- ILLER 

ENTER!-  RESIN  RESISTIVITY 

♦  UE’M6j 

ENTER!-  FIBER  RESISTIVITIES 

( LONGITUDINAL  .. TRANSVERSE  ) 

♦  7  E  ~  6  >  1  ♦  7'  c!  “  6  ? 

ENTER:>  FIBER  VOLUME 
START  VALUE? 

4  0»  1.0?  0.>02? 

enter;::-  aspect  RATIO  OF  FIBER 

0.0? 

ENTER 


ENTEF’-:;>  HERMANS  OR  I  EiM  TAT  I  ON  PARAMETER 


TYPE  OF  AVERAGING 

( 1 >  IN-PLANF  ?  TWO-D I  MENS IDNAL 

(2)  A.KIAL  SYMMETRIC?  THREE  — DIMENS IONA 


FRACTIONS  TO  SCAN 
END  VAI..UE?  INCREMENT 


.  >-■  ? 


143 


Table  9  (continued) 


DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<DATA  CHECIO 
CALIBRATION  PARAMETERS 
ALPHA  =  0.4100 

,  beta  =  0.2900 


}•  DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 
H 


•f  DO  YOU  WISH  A  TABLE  OF  DATA? 

Y 

+  DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 

ENTER>  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 

(1)  Kll 

( 2 )  K22 

(3)  K33 

1  9 


Table  9  (continued) 
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ELECTRICAL 

RESITIv'ITIES  OF  COMPOSITE  MATERIAL  I  AXIAL  CASE 


'•/F 

F 

Kll 

K22 

!\  -.i 

0,0200 

0,5000 

4,2314E+15 

1 ,3083E-M6 

1 ,3083E+16 

0 , 0400 

0 ♦ 5000 

2, 1348E+15 

7,3967Efl5 

7,3967E+15 

0,0600 

\j .  5  0  0  0 

1 ,3883E-rl5 

5.0266Efl5 

5,0266EF15 

0  +  0800 

0,5000 

1 ,0062ET15 

3 ,7291Efl5 

3,7291E+15 

0.1000 

0 , 5000 

1 ,5154E-l-07 

6,061  SET 07 

6,0618E+07 

0,1200 

0,5000 

2,1250E-05 

8, 500 OF -05 

8 ,5000E-05 

0, 1400 

0 , 5000 

1 ,8214E-05 

7,2357E-05 

7. 2357 E- 05 

0 . 1600 

0,5000 

1 ,5933E-05 

6,3750E-05 

6,3750E-05 

0,1800 

0 . 5000 

1 ,4167E~05 

5,6667E-05 

5,6667E-05 

0,2000 

0 , 5000 

1 ,2750E-05 

5, 1000 E- 05 

5 . lOOOE-05 

0.2200 

0,5000 

1  , 1591E-05 

4,6364E-05 

4, 636 4E -05 

0,2400 

0.5000 

1 .0625E-05 

4,2500E-05 

4,2500E-05 

0.2600 

0 , 5000 

9,8077E-06 

3.9231E-05 

3 .9231E-05 

0,2300 

0,5000 

9,1071E-06 

3,6429E-05 

3,6429E-05 

0,3000 

0,5000 

8 . 5000E“-06 

3 . 4000E-05 

3,4000E-05 

0,3200 

0,5000 

7.9688E‘-()6 

3,lS75E-05 

3, 1S75E-05 

0  ♦  *./ 

0.5000 

7,5000E-06 

3,OOOOE-05 

3  .  OOOOE-05 

0 . 3600 

0,5000 

7,0333E-06 

2,8333E-05 

2,8333E-05 

0.3800 

0.5000 

6,7105E--06 

2,6842E-05 

2 .6342E-05 

0.4000 

0,5000 

6.3  7 5  0  E,  ““  0  6 

2,:5500E-05 

2.5500E~05 

0 . 4200 

0 . 5000 

6,0714E-06 

2,4206E-05 

2.42S6E-05 

0,4400 

0 , 5000 

5.7935E-06 

2,3132E-05 

2,3182E-05 

0 . 4600 

0,5000 

5,5435E-06 

2 ,2174E-05 

2 . 2 1  7  4  E  0  5 

0 , 4800 

0,5000 

5,3125E-06 

2,1250E-05 

2,1250E-05 

0 , 5000 

0 , 5000 

5,1000E-06 

2,0400E-05 

2.0400E~05 

0,5200 

0 , 5000 

4,903BE-06 

l,9615E-05 

1 ,9615E-05 

0 , 5400 

0 , 5000 

4,7222£~06 

1 ,8889E~05 

1 ,833 9E -05 

0,5600 

0 , 5000 

4,5536E--06 

1 .8214E-05 

1 ,821 4E -05 

0,5800 

0,5000 

4 ,3966E-06 

1 .7586E-05 

1 ,7586E~05 

0 , 6000 

0,5000 

4,2500E-06 

1 ,700011-05 

1 ,7000E-05 

0 , 6200 

0 , 5000 

4, 1129E-06 

1 ,6452E~05 

1 ,6452E-05 

0 , 6400 

0 , 5000 

3,9844E-06 

1 .5938E-05 

1 ,5938E-05 

0,6600 

0 , 5000 

3,8636E-06 

1  ,5455E-05 

1 .5455E-05 

0,6800 

0,5000 

3 , 7500E-06 

1  ,5000E-05 

1 ,5000E-05 

0,7000 

0,5000 

3,6429E-06 

1 ,457 IE- 05 

1 ,4571E-05 

0,7200 

0 , 5000 

3,5417E-06 

1 ,4167E-05 

1 .4167E-05 

0,7400 

0,5000 

3,4459E-06 

1 ,3784E-05 

1 ,3734E-05 

0 ,7600 

0 , 5000 

3,3553E--06 

1 ,3421E-05 

1 ,3421E-05 

0,7800 

0,5000 

3, 2692 E -06 

1 ,3077E-05 

1 .3077E-05 

0 , 8000 

0 , 5000 

3, 1 37 5 E -06 

l,2750E-05 

1 ,2750E-05 

0,8200 

0,5000 

3 , 1098E-0c> 

l,2439E-05 

1 ,2439E-05 

0 ♦ 3400 

0,5000 

l,2143E-05 

1 ,2143E-05 

0,8600 

0,5000 

2,9651E-06 

1 , 1860E-05 

1  ,  1860E-05 

0 ,8800 

0,5000 

2,8977E-06 

1 , 1591E-05 

1  ,  1591E-05 

0 , 9000 

0 , 5000 

2 , 8333E-06 

1 ,1333E-05 

1  ,  1 333E-05 

0 , 9200 

0 , 5000 

2,7717E-06 

1 , 1087E-05 

1 , 1087E-05 

0,9400 

0,5000 

2,7128E-06 

1 ,0851£-05 

1 .085  IE -05 

0 , 96'00 

0 , 5000 

2,6563E-06 

1 ,0625E-05 

1 , 0625E-05 

0,9800 

0,5000 

2,6020E-06 

1 ,0403E-05 

1 ,0403E-05 

1 , 0000 

0 , 5000 

1 , 7 4  6  y  E — 0  6 

1 ,3223E-06 

1 .8223E-06 

Tcible  9  (continued) 
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1  READY  FOR  PLOT? 
Y 


ENTER>  VERTICLE  AXIS  HODE 
1  (1)  PLOT  LINEAR  SCALE  FOR  LOGlO(Kii) 

(2)  PLOT  LOG  CYCLE  SCALE  FOR  LOGIC (K:i;i) 

1  S' 

ENTER>  GRAPHICS  OUTPUT  DEUICE  ■  ' 

(1)  TERMINAL 

(2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER  ‘ 

<4)  END  PROGRAM 


ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  yAI..UEf  END  VALUE?  INCREMENT 

0,0f  I4O?  0*02? 

ENTER>  ASPECT  RATIO  OF  FIBEI'^S 
50*0? 

ENTER>  TYPE  OF, AVERAGING 

(1)  IN-PLANE  5  TUO-DIMENSIONAL 

(2)  AXIAI.  SYMMETRIC?  THREE-DIMENSIONAL 


ENTER>  HERMANS  ORIENTATION  PARAMETER 
0*5? 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<DATA  CHECIO 
CALIBRATION  PARAMETERS 
ALPHA  =  0*4100 

BETA  =  0*2900 


■f  DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 
N 

■f  DO  YOU  WISH  A  TABLE  OF  DATA? 

N 

•f  DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 


Table  9  (continued) 
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ENTER>  GF-:APHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 
(1)  Kll 
■;2)  K22 
(3)  K33 


}•  READY  FOR  PLOT? 

Y 

ENTER>  'v'ERTICLE  AXIS  MODE 

(1)  PLOT  LINEAR  SCAi  E  FOR  LOGlOCKii) 

■  (2)  PLOT  LOG  CYCLE  SCALE  FDR  LOGIC  (K;l 
;i.  i- 

GRAPHICS  OUTPUT  DEVICE 
<1)  TERMINAL 

(2)  FLATBED  PLOTTER 


REPEAT  OPTION 

(1)  REPEAT  SAME  PL.OT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 
<3>  START  OVER 

(4)  END  PROGRAM 


ENTE:R>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE  f  END  VALUE  i.  INCREMENT 

0,0?  1,0?  0.02? 

ENTER?  ASPECT  RATIO  OF  FIBERS 
50,0? 

ENTER?  TYPE  OF  AVERAGING 

(1)  IN-PLANE?  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC?  THREE-DIMENSIONAI 

2? 

ENTER?  HERMANS  ORIENTATION  PARAMETER 
0,5? 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

CDATA  CHECK? 

CALIBRATION  PARAMETERS 
ALPHA  -  0,4100 

BETA  =  0,2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 


i  DO  YOU  WISH  A  TABLE  OF  DATA? 
N 

■f  DO  YOU  WISH  A  COMPUTER  PLOT? 


ENTER? 


ENTER? 


Table  9  (continued) 
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ENTER>  GRAF!-! I C  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 

(1)  Kll 

(2)  K22 

(3)  K33 
3  ? 

■f  READY  FOR  PLOT? 

Y 

ENTER>  VERTICLE  AXIS  HODE 

(1>  PLOT  LINEAR  SCALE  FOR  LOGlO(Kii) 
<2>  PLOT  LOG  CYCLE  SCALE  FOR  L0G10<Ki 
;l  1- 

ENTER>  GRAPHICS  OUTPUT  DEVICE 

(1)  TERHINAL 
<2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAHF.  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER 

(4)  END  PROGRAM 
4  > 


Figure  8.  Electrical  resistivity  for  copper  fibers  and  PPS  resin,  transverse  direction. 
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7)  Copper  Flakes  and  PPS  Resin 


Copper  Flakes^ 

Resin 

Resistivity  (ohm-cm) 

1.70x10“^ 

4.5x1016 

Volume  Fraction 

Variable 

Aspect  Ratio 

1 .OxlO"5 

— 

Orientation  Parameter , 
fa  ,  (axial) 

1.0 

Executions  are  given  in  Table  10  for  the  computer 
graphs  depicted  in  Figures  10,  11  and  12  for  ,  K2  , 
and  K3  ,  respectively. 


Note;  The  influence  of  geometrical  irregularities  and/or 
surface  coatings  may  be  taken  into  account  for 
replacing  the  filler  resistivity  by  the  effective 
resistivity. 

^Effective  resistivity  of  copper  =  1.70x10-3  ohm-cm. 
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Table  10* 

Calculations  for  Resistivity  for 
Copser  Flakes  and  F’F'S  F'lesin* 


CENTER  FOR  COMPOSITE  MATERIALS 
COMPUTER  SOFTWARE  SYSTEMS 

H I C R 0 ME C H A H 1 C A L  M 0 D EL  F 0 1^  C A L C U L A T 1 N G  i 
COMPOSITE  MATERIAL  E1..ECTICAI,.  AND  THERMAL 
RESISTIUITIES  AND  CONDUCTIOITES 

UERSIONi  051 1S2 


ENT'EFD-  CALCULATION  OPTION 

<1)  CAi..CULATE  EL.ECTRICA!..  PROPERTIES 
(2)  CAi  CULATE  THFRMAI.  PROPERTIES 


ENTER>  INPUT  OPTION 

(1)  INPUT  RESISTIVITY  DATA 

(2)  INPUT  CONDUCTIVITY  DATA 


f  DO  YOU  WISH  TO  INCLUDE  A  FILLER  PHASE? 

ri 

TWO  PHASE  SYSTEM  ASSUMED J  RESIN/F ILLER 

ENTER>  RESIN  RESISTIVITY 
4,5Etl6* 

ENTER>  FIBER  RESISTIVITIES 

(L0N6ITUDINAI..?  TRANSVERSE) 
l,7E-6f  1.7E-6? 

ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE  END  VALUE?  INCREMENT 
v+0>l*0?0*02> 

ENTER>  ASPECT  RATIO  OF  FIBERS 
1 ♦OE-S? 

ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE»  TWO-DIMENSIONAL 

<2)  AXIAL  SYMMETRIC)  THREE-DIMENSIONAL 


ENTER>  HERMANS  ORIENTATION  PARAMETER 

1,0? 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

CD  AT  A  check;:-- 
CALI  BRAT I  ON  PARAMETERS 
AI..PHA  =  0,4100 

BETA  =  0,2900 
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Table  10  (continued) 


t  DO  YOU  WISH  TO  CHANGE  THESE  UALUl 
N 

t  DO  YOU  WISH  A  TABLE  OF  DATA? 

Y 

•f  DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 

ENTER>  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 
<1>  Kll 

( 2 )  K2? 

(3)  K33 

1? 


Table  10  (continued) 


ELECTRICAL 


lESITIVITIES 

OF  COHF'OS 

ITE  MATERIAL  1 

AXIAL  CASE 

VF 

F 

Kll 

K22 

K33 

0,0200 

1 , 0000 

4,4100E-M6 

1 .6965E-H3 

^  5  F  -f  1 3 

0 , 0400 

1,0000 

4,3200E+16 

e,1414Ef 12 

3+ 1414Efl2 

0,0600 

1,0000 

4.2300E+16 

5 ,2042i:  +  12 

5.2042E+12 

0,0800 

1.0000 

4, 1400E+16 

3,7389Efl2 

3 .7389E+12 

0,1000 

1,0000 

4  ♦  OSOOE'flo 

2.3626Efl2 

2.S626E+12 

0,1200 

1,0000 

3.9600E+16 

2,2S07E+12 

2.2807E+12 

0,1400 

1,0000  , 

3,8700ET16 

1,8670E-H2 

1,8670E+12 

0,1600 

1,0000 

3.7800E+16 

1 .5586E+12 

1  +  5  5  3  6  E  f  1 2 

0,1800 

1 , 0000 

3,6900E+16 

1 ,3202E+12 

1 +3202Eil2 

0,2000 

1,0000 

3 ♦ 6000E+ 1 6 

1 , 1309E-)-12 

1 , 1309E+12 

0*2200 

1.0000 

3,5100E-M6 

9,7736EH1 

9,7736E+11 

0 , 2400 

1  ♦  0000 

3,4200ET16 

8,5056E-f-ll 

8,5056£+ll 

0+2600 

1 ♦ 0000 

3,3300EH6 

7 ,4436£  fll 

7,4436E+11 

0 , 2800 

1 , 0000 

3,2400E+16 

6 , 5433Ei 1 1 

6  *  5433E  f 1 1 

0,3000 

1.0000 

3,1500Efl.6 

5,772511  +  11 

5+7725Ef 11 

0,3200 

1 ,0000 

3,0600E"}-16 

5,1069E+11 

5, 1069E+11 

0 ♦ 3400 

1 ♦ 0000 

2,9700E+16 

4.5280E-M1 

4,5280E+11 

0  +  3600 

1 .0000 

2,8800E-!-16 

4,0212E+11 

4.0212E+11 

0,3800 

1 , 0000 

2,7900E+16 

3,574SE+11 

3.5748E+1 1 

0,4000 

1 , 0000 

2,7000E-+16 

3. 1723E  +  11 

3,1723E+11 

0,4200 

1  +  0000 

2,6100Efl6 

2,7323E+11 

2,7323E+1 1 

0,4400 

1 , 0000 

2.5200E+16 

2,0181E+1 1 

2.0181E+11 

0 , 4600 

1.0000 

2.4300Efl6 

1 ,0102E+1 1 

1 ,0102E+11 

0,4800 

1 ,0000 

2,3400E-H6 

2.7333E+10 

2 , 7333E  +  10 

0,5000 

1 , 0000 

2,2500E+16 

3,4660E+09 

3  ♦  4-660E+09 

0,5200 

1 , 0000 

2,1 600ET 1 6 

1 .90S6E+08 

1 ,9086E+08 

0,5400 

1 . 0000 

2.0700Et16 

4,3833E+06 

4.3833Ef06 

0 , 5600 

1 , 0000 

1 ,9800E+16 

4, 1042E+04 

4 , 1042E+04 

0,5800 

1 , 0000 

1 ,8900E+16 

1  ,5485E+02 

1 ,5485E+02 

0 , 6000 

1 , 0000 

1 ,8000E+16 

2 , 8  -6  o  iz. 0 

2 ,3333E-06 

0,6200 

1 , 0000 

1 ,7100E+16 

2,7419£-0.6 

2,7419E-06 

0,6400 

1,0000 

1 ,6200E+16 

2,6563E-06 

2,6563E-06 

0 , 6'600 

1,0000 

1 ,5300E+16 

2,5758E-06 

2.5758E-06 

0 , 6-800 

1 , 0000 

1  ,4400E+16 

2,5000E-06 

2.5000E-06 

0,7000 

1 , 0000 

1 ,3500E-fl6 

2,42Ei6E-0-6 

2,4286E-06 

0 , 7200  ■' 

1,0000 

1  .2599E  +  16 

2,361 lE-06 

2 , 3  -6 1 1 E  -  0  6 

0,7400 

1,0000 

1,1699E-H6 

2,2973E-06 

2,2973E-06 

0 ,7-600 

1,0000 

1 .0799E-fl6 

2,2368E-06 

2,236SE-06 

0 , 7800 

1,0000 

9,3994E-1-15 

2. 1795E-06 

2,1795E~06 

0,8000 

1.0000 

8,9994E-tl5 

2 . 1250E-06 

2, 1250E-06 

0,8200 

1,0000 

8,0994ET15 

2,0732E-06 

2 .0732E-06 

0,8400 

1,0000 

7,1994Efl5 

2,0238E-06 

2,023SE-06 

0 , 8600 

1 , 0000 

-6,2994E+15 

1 .9767E-06 

1 .9767E-06 

0,8800 

1 , 0000 

5,3994E+15 

1 ,9318E-06 

1 ,9318E-06 

0 , 9000 

1,0000 

4,4994Efl5 

1  ,3S39E-06 

1 .8839E-06 

0,9200 

1,0000 

3,5993E+15 

1 ,8478E-06 

1 ,S478E-06 

0,9400 

1 . 0000 

2,6993E-!-15 

1 , 8  0  8  o  E 0 

1 ,S035E-06 

0 , 9600 

1,0000 

1 ,7993E-}-15 

1 ,7708E-0-6 

l,7703E-06 

0,9800 

1,0000 

8,9931E+14 

1 .7347E-06 

1 ,7347E-06 

1,0000 

1 , 0000 

9.4805E-03 

1 , 7000E-06 

1 , 7  0  0  0  £  -  0  -6 
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•f  READY  FOR  PLOT? 

Y  ■  • 

ENTER>  VERTICLE  AXIS  HODE 

(1)  PLOT  LINEAR  SCALE  FOR  LOGIC (Kii) 

<2)  PLOT  LOG  CYCLE  SCAI-E  FOR  L0G10<Kii) 

If 

ENTER>  GRAPHICS  OUTPUT  DEVICE 
(1>  TERMINAL 

(2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEU  MATERIAL  GEOMETRY 
<3)  START  OVER 

<4)  END  PROGRAM 


ENTER>  FIBER  VOLUME  FRACTIONS  TO  SCAN 

START  VALUE »  END  VALUE?  INCREMENT 

0»0?  1.0?  0,02? 

ENTER>  ASPECT  RATIO  OF  FIBERS 
;l.  .OE-5? 

ENTER>  TYPE  OF  AVERAGING 

(1)  IN-PLANE  5  TWO-DIMENSIONAL 

(2)  AXIAL  SYMMETRIC*?  THREE-DIMENSIONAL 
2?  ' 

ENTER>  HERMANS  ORIENTATION  PARAMETER 

1.0? 

DEFLAULT  ELECTRICAL  CALIBRATION  PARAMETERS 

<DATA  CHECK> 

CALIBRATION  PARAMETERS 
ALPHA  =  0.4100 

BETA  =  0.2900 


DO  YOU  WISH  TO  CHANGE  THESE  VALUES? 
N 

DO  YOU  WISH  A  TABLE  OF  DATA? 

N 

DO  YOU  WISH  A  COMPUTER  PLOT? 

Y 

ENTER>  GRAPHIC  OUTPUT  DATA 

SELECT  COORDINATE  DIRECTION 
(1)  Kll 
(2>  K22 

(3)  K33 
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ENTER>  GRAPHIC  OUTPUT  DATA 


SELECT  COORDINATE  DIRECTION 

(1) 

Kll 

(2) 

K22 

(3> 

K33 

3  f 

READY  FOR  PLOT? 

Y 

ENTER>  VERT I CL E  AXIS  HODE 

(1)  PLOT  LINEAR  SCALE  FOR  LOGKXKii) 

(2)  PLOT  LOG  CYCLE  SCALE  FOR  LOGlOCKi 
If 

ENTER>  GRAPHICS  OUTPUT  DEVICE 
(1)  TERMINAI.. 

•:2)  FLATBED  PLOTTER 


ENTER>  REPEAT  OPTION 

(1)  REPEAT  SAME  PLOT 

(2)  ENTER  NEW  MATERIAL  GEOMETRY 

(3)  START  OVER 
<4)  END  PROGRAM 

4f 


Figure  10.  Electrical  resistivity  for  copper  flakes  and  PPS  resin,  longitudinal  direction. 
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Electrical  resistivity  for  copper  flakes  and  PPS  resin,  transverse  direction. 
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Figure  12.  Electrical  resistivity  for  copper  flakes  and  PPS  resin,  perpendicular  direction. 
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